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Within the framework of the hard thermal loop effective theory we derive a system of 
Boltzmann-like kinetic equations taking into account the simplest processes of nonlinear 
interaction of soft fermionic and bosonic QCD plasma excitations: elastic scattering of 
soft- (anti) quark excitations off soft-gluon and soft-quark excitations, pair production of 
soft quark-antiquark excitations, annihilation into two soft-gluon excitations. The matrix 
elements of these processes to leading order in the coupling constant g are obtained. The 
iterative method of calculation of the matrix elements for the higher processes of soft- 
mode interactions is proposed. The most general expression for the emitted radiant power 
induced by the effective currents and effective sources in a quark-gluon plasma (QGP) 
taking into account an existence of fermion sector of plasma excitations is defined. The 
explicit form of the linearized Boltzmann equation accounting for scattering of color(less) 
plasminos off color (less) plasmons is written out. 
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1 Introduction 



In our previous papers PJ E] on the basis of a purely gauge sector of the Blaizot-Iancu 
equations |3] complemented by the Wong equation jlj we have studied in great detail 
the nonlinear dynamics of soft boson excitations in hot QCD plasma. In this and two 
accompanying papers we would like to enlarge the analysis carried out in Refs. P3I2] on the 
fermion sector of soft plasma excitations. For solving this problem it is necessary to use 
the whole system of the Blaizot-Iancu equations. It is well known that hot QCD plasma 
including massless quarks (and antiquarks) possesses doubled dispersion relation for soft 
fermion plasma excitations [3]. The first branch describes normal particle excitations (in 
the subsequent discussion designated by a symbol "+") with relation between chirality and 
helicity at zero temperature. The second branch is purely collective excitation (designated 
further by symbol " — "), where usual relation between chirality and helicity is flipped 
jH] • It has been called a plasmino [7j by analogy with the plasmon (longitudinal) mode of 
gluons having also a purely collective character. 

Here, in the first part we have restricted ourselves to study of the processes connected 
with interaction of the soft fermion and boson excitations among themselves without ex- 
change of energy with hard thermal (or external) partons. We consider the processes 
of nonlinear interaction of soft excitations within the framework of the kinetic approach 
by the determination of corresponding set of coupled Boltzmann-like kinetic equations 
describing some non-linear aspects of the dynamics of the soft collective excitations in 
high-temperature QGP. The construction of such equations is generally straightforward 
at a conceptual level, but it is tedious and intricate in practice. By virtue of complex- 
ity and awkwardness of the kinetic equations for soft plasma modes taking into account 
simultaneously the fermion and boson degrees of freedom, we will consider in details 
only the simplest processes of the nonlinear interaction: soft- (anti) quark -soft-gluon and 
soft- (anti) quark -soft-quark elastic scattering, production of soft quark-antiquark pair by 
fusion of two soft gluons, and annihilation of soft quark-antiquark pair into two soft-gluon 
excitations. At least for weakly excited system corresponding to the level of thermal fluc- 
tuations at the soft momentum scale these processes are the basic those of the nonlinear 
interaction of soft excitations in the medium. In deriving the kinetic equations for soft 
fermion and boson modes we use oversimplified approach. Proceeding from the semiclas- 
sical character of the problem, we initially assume that the structure of the collision terms 
is determined by Fermi's golden rules and thus we focus our efforts exclusively on cal- 
culating the scattering probabilities. We calculate these probabilities in a direct fashion 
by simple extracting all relevant contributions in scattering processes 2 — > 2 enumerated 
above. We don't give an explicit proof of a gauge invariance of the scattering probabilities 
obtained. This can be performed in spirit of our previous paper [H] by using the effective 
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Ward identities for hot gauge theories HH] ■ 

It should be particularly noted that isolation and study of the nonlinear dynamics of 
soft-mode interactions only in QGP presented in this work have sufficiently artificial char- 
acter and bear rather limited sense. Really, by the efforts of a number of authors (see, for 
example |TTJ^]) it was shown that the dominant interactions in QGP are generally those 
between the soft modes of interest and the hard particles in the thermal bath. Moreover, 
the general transport or equilibration properties of the plasma are in fact controlled by the 
hard particles (those with typical momenta of order T, the temperature of the system). In 
our following papers the interaction processes of soft and hard modes of different statistics 
among themselves and also radiation processes of soft modes induced by collisions of hard 
thermal particles will be considered in details. The construction of the effective theory for 
such processes has shown that it is more convenient to consider separately a theory for 
simpler (in the conceptual plan) interaction processes of the soft modes. Matrix elements 
of interaction processes of soft and hard excitations include, as constitutive blocks, the 
expressions defining matrix elements of self-interaction of soft excitations provided that 
external momenta are put on relevant mass-shells. This significantly facilitates general 
analysis and classification of individual terms in very complicated and tangled expressions 
for matrix elements of the interaction processes between the soft excitations and the hard 
particles including simultaneously the soft and hard modes of Fermi and Bose statistics. 

Incidentally there are circumstances at which the kinetic equations for the soft quasi- 
particles may be relevant independently: for instance, when we are interested in the 
thermalization of such soft modes or in their damping rates. Besides at studying the 
processes with soft-gluon excitations only we have shown PJ E| that there exists a certain 
level of intensity of plasma excitations above which the nonlinear interaction processes 
of soft boson excitations among themselves start to dominate over the interaction pro- 
cesses of soft modes with hard thermal particles, thus defining in particular qualitatively 
different mechanism for the damping of a soft gluons. In other words, for large values 
of soft-gluon occupation number (i.e. when a subsystem of soft boson excitations is far 
from equilibrium) dynamics of hot QCD plasma is defined by scattering processes of soft 
quasiparticles off each other. In the second part [T3] we will show that and for the case of 
fermion sector of plasma excitations in QGP there exists a critical threshold above which 
the interaction processes of soft fermion and boson excitations among themselves start to 
play dominating role in general dynamics of the system. 

In the light of discussed above a study of soft-excitations dynamics only can represent 
independent interest, however this is beyond the scope of our work. Our ultimate goal is 
the study of new (not considered earlier in the literature) possible mechanisms of energy 
losses of high-energy particles propagating through hot QCD medium associated with an 
existence in medium soft- (anti) quark excitations on a level with soft-gluon excitations. 
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This is a main subject of the following papers. Research of the energy losses of energetic 
partons in QGP at present is of great interest with respect to jet quenching phenomenon 
[T4*l IT5] . The basic purpose of present work is development of the general technique for 
calculation of relevant probabilities and in particular, deriving the semiclassical formula 
for the emitted radiant power in QGP taking into account fermion sector of plasma 
excitations (section 8). This formula (with the further updating) is basic formula for 
defining energy losses of energetic parton (quark or gluon) travelling through quark-gluon 
plasma. 

The derivation of the Boltzmann equation describing evolution of the number density 
of fermion quasiparticles in QGP was also considered by Niegawa JH] • Within the frame- 
work of the Keldysh-Schwinger formalism the kinetic equations for normal and abnormal 
modes here, are defined from the requirement of the absence of the large contributions 
due to pinch singularities of the perturbative scheme proposed in Ref. [lfij . The formalism 
suggested by Niegawa is quite rigorous and possesses great generality at least for quasi- 
uniform systems near equilibrium. Unfortunately in this paper particular expressions for 
collision terms wasn't given. In this sense, our work having more phenomenological char- 
acter supplements the paper jTHj with examples of concrete calculation of the collision 
terms within the HTL effective theory. 

The paper is organized as follows. In section 2, preliminary comments with regard 
to deriving the system of the Boltzmann equations describing the processes of nonlinear 
interaction of soft fermion and boson excitations in QGP are given. In section 3, the 
system of nonlinear integral equations for gauge potential A ^ and quark wave function 
ip l a is written out. On the basis of perturbative solutions of these equations the notions 
of the effective currents and sources playing a key position in our research are given. In 
section 4, an algorithm of the successive calculation of the effective amplitudes defining 
matrix elements of the processes we are interested is presented. In sections 5, 6 and 7, the 
details of calculations of the probabilities for the simplest processes of interaction: soft- 
(anti) quark -soft-gluon and soft- ( ant i) quark -soft- quark elastic scattering, production of 
soft quark-antiquark pair by fusion of two soft gluons, and annihilation of soft quark- 
antiquark pair into two soft-gluon excitations are presented. Section 8 is devoted to 
deriving the most general formula for the emitted radiant power in QGP. It was shown 
how this formula within the framework of Tsytovich correspondence principle allows by a 
direct way to calculate the probabilities of higher order processes of nonlinear interaction 
of soft plasma modes. In section 9, the explicit form of linearized Boltzmann equation 
for colorless plasminos is written out and some speculations concerning the extension 
of the formalism developed to colorcharged plasminos are given. In Conclusion some 
features of dynamics of soft excitations is briefly discussed and connection with more 
traditional approaches based on imaginary- and real-time field theories is traced. Finally, 
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in Appendixes summary of various properties of HTL-amplitudes used in the text and 
the explicit form of the conjugate effective amplitudes are given. 



2 Preliminaries 



In this paper we consider a construction of effective kinetic theory for colorless soft- 
quark and soft-gluon plasma excitations propagating in a hot quark-gluon plasma. We 
assume that localized number densities of soft-quark n ± (q, a;) = (n* y ), soft-antiquark 
fr^q, x) 
space 



n^i), and soft-gluon 7V*''(k, x) = (N^ lab ) excitations are diagonal in a color 



n 



±ij 



5 tj n± , ij = 5 tj n± , iV£ 1 ab = 5 ab N^ 1 , 



where i, j = 1,2, . . . , N c and a, b = 1, ... , N 2 — 1 for SU (N c ) gauge group. We consider a 
change of the number densities of the colorless soft- (anti) quark excitations as a result of 
their interactions among themselves and with the colorless soft-gluon excitations. 

The dispersion relations for soft-quark modes c<j ± (q) = Wq and soft-gluon modes 
uj t,l (k) = lu^ are defined by equations 



Re*A±V,q) = °> Re*A- 1M (A; ,k) = 0, 



(2.1) 
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are inverse quark and gluon scalar propagators. cJq = g 2 CpT 2 /8 and o;^ = g 2 {2N c + 
Nf)T 2 / 18 are plasma frequencies squared of the quark and gluon sectors of plasma excita- 
tions. Hereafter, we denote momenta of the soft-quark fields by q, q', q±, . . . and momenta 
of the soft-gauge fields by k, k', k 1: . . . . 

We expect a time-space evolution of the scalar functions f — ± and b = t, I 
to be described by the self-consistent system of Boltzmann-like equations 



^? +v (/) 
dt q 



dNj b) 



dt 



<9x 
<9x 



J- = -< } r d W K^1 + (1 - n^KN^ 1 ], 



(2.2) 
(2.3) 



where = du^/dq and v k = du^ /dk are the group velocities of soft fermionic and 
bosonic excitations, respectively. The generalized decay rates T^/' ^ and inverse decay rates 
p(/>&) are ( non _ii near ) functionals of the soft- (anti) quark and soft-gluon number densities. 
Here, for the sake of brevity we drop dependence on soft-antiquark occupation numbers 
on the right-hand side of Eqs. ()2.2)1 and (|2.Hj) . The equation for is obtained from 



n 



(|2.2j) by replacement 



The decay and regenerating rates can be formally represented in the form of functional 
expansion in powers of the soft- (anti) quark and soft-gluon number densities 

oo oo 

rf VX'*] = Er? 6)(M) [nJX''], r^[n^^ 1 ) = ErF' 6)(2 " + VX1, 

n=l n=l 

(2-4) 

where T < "/'^ 2n+1 ' > [n^ , N^ 1 ] collect the contributions of the total nth power in n^, and 
N± 1 . Unfortunately, the general structure of the expressions for arbitrary n is very cum- 
bersome and therefore we restrict our consideration to the simplest case for n = 1. 

The fermion decay and regenerating rates in the lowest order of the nonlinear inter- 
action (n = 1 in Eq. ()2.4j) ) can be formally represented in the following form: 



rJ/VjX-'HE E 2/<^/ l62) ^ l62) (q,k i; qi,k 2 )(l-n^X^(l+iV k2 



/l=± b!,b2=t,l 



(62) 

(2.5) 



1" / Ulqg^gg UJgg^gg 



Q3 I I ' 



dquki.^fl + JV^Jl + ^J 

+ E { / </iV/ 2/3) <^ /2/3) (q, Qi; q 2 , <b)n™ (1 - < 2) )) (1 - nj{) 

h,h,h=± < J 

+ 2jdJ q ^ f3) w^ f3 \q, qi ; q 2j q 3 )n^(l-^)(l-^) 

and in turn, 

rFVt =E , E j2/<^/ l62) <^ /l62) (q,k i; qi ,k 2 )<;)(l+<))iv/^ 

+ E I jdr^gf^w^^^^ q 2 ,q 3 )(l 

+ Z / alqg^qq Wqq^qq Ql, q2, 



/l=± b!,b 2 =t,l 
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4- UT {ffl ' blb2) in Uh-Mb2) l 

T" / Ulgq^gg w qq ^g g 



r) (/l) N \ r) (/2) r ,(/3) 
"qi ,Pq2 ' 4 q3 



qi / q2 q3 



The functions w^ b ^' q { / lb2 \ Wg^' g g lb2 ^ are probabilities for 'elastic' scattering of soft-quark 
excitations off soft-gluon excitations and annihilation of soft quark- ant iquark pair into two 
soft-gluon excitations, and w^' q ^ 3 \ Wqq^qq 2 ^ are probabilities for 'elastic' scattering 
of soft-quark excitations off soft-quark and soft-antiquark excitations, respectively. The 
phase-space measures for these processes are 

/ <W> - /<§3 /(^3 j^i (S-rMq+k.-q, - fc) -«£>), 

(2.7) 

/^■«. / *l | / *l (2,)^( q+qi -k 1 -k 2 ) 5 (^) + .<«_4M_4« ) , 

= S /*L/*L/^ + * - q2 - *) 

x^+ w («- w 2 2) -#). 

For boson sector of the plasma excitations the generalized rates and rf ^ to the 
lowest order in the nonlinear interactions of soft modes have a similar structure 



r?KX''l = E E /^ ( S l/2) « /2) (k, qi ; k 1)q2 )(i+0^)(i-n^)) 

&1=M /l,/2=± l J 

(2.8) 

+J<n£*£to *<£#*>(k, q i; kl , q 2 ) (l + (l - «W 



n q 2 



rf Vt ^'']=E E { / ^S l/2) « /2) ( k . oi; k i, q 2 )< l) (i - 

&!=*,/ /l,/2=± ^ 

where, w^JX'gq^\ w gq-^gq^ are probabilities for 'elastic' scattering of soft-gluon excita- 
tions off soft-quark and soft-antiquark excitations, respectively, and Wg^Jq^ is the prob- 
ability for soft quark-antiquark pair creation by two soft-gluon fusion. The phase-space 
measures for these processes are 

/«« = * /|5L/^ (2 ^ (k + qi _ kl _ , 2) (, 9) 
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xirf+<''-<' ) -^'), 

We have written out above the most general expressions for the fermion and boson 
decay and regenerating rates to the lowest order in the nonlinear interaction taking into 
account all possible channels of transitions from initial two-quasiparticle states to different 
types of the final two-quasiparticle states. However it is clear that not all these processes 
of scattering, annihilations and fusions are kinematically permissible 1 . In this work we 
do not set as a purpose to determine all permissible channels of reactions that require as 
a first step a study of a system of equations determining the conservation laws of energy 
and momentum of two-quasiparticles interacting system. Unfortunately at present there 
not exist general analytic methods 2 permitting in a direct way to define whether or not 
^-functions have supports different from zero in integration measures ()2.7j) and ()2.9|) and 
thus for solving this problem it should be used numerical methods. In the remaining 
part of the paper we restrict our consideration to calculation of the probabilities of only 
processes, which certainly exist (for example, the elastic scattering) leaving the general 
case for a separate research. 



3 Effective currents and effective amplitudes 

We use the metric g^ u = diag(l, — 1, — 1, — 1), choose units such that c = k B = 1 and 
note x = (t, x), k = (fc ,k), q = ((fo,q) etc. As was mentioned above, we consider 
SU(iV c ) gauge theory with Nf flavors of massless quarks. The color indices for the adjoint 
representation a,b, . . . run from 1 to N^ — l, while those for the fundamental representation 
. . . run from 1 to iV c . The Greek indices a,/3, . . . for the spinor representation run 
from 1 to 4. In the following discussion we will use notation adopted in our previous 
paper jH]. 

The input equations for the effective theory under consideration are self-consistent 
system of the field equations: the Yang-Mills equation for gauge potential (k) (Eq. (4.9) 

^^Here, as a kinematic permissibility we mean generally speaking not only the existence of a solution 
for the system of conservation laws of energy and momentum corresponding to concrete reaction of the 
nonlinear interaction of two quasiparticles, but also selection over helicity conservation of a quasiparticle 
system. 

2 In the work of Kadomtsev and Kontorovich ^7] the geometrical method of solving the system of 
conservation laws for two-to-two scattering processes of quasiparticles with arbitrary dispersion laws was 
suggested. Unfortunately, this approach having sufficiently general character is restricted only to the case 
of two-dimensional momenta of quasiparticles. 
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in Ref. 0) 

*V^(k)A™(k) = -j^ 2)a (AA)(k) - j^ a (A,A,A)(k) -f^ a $^){k) 

-f^\A^){k), (3.1) 

and the Dirac equation for soft-quark field (Eq. (3.3) in jH]) supplemented by its Dirac 
conjugate equation for completeness of the description 

*S^(q)4(q) = -^ 1)l (A^)(q) -^(AA^iq), (3.2) 

*S£(-q) = $){-q) + fj^A*, A*,j)(-q). 

Here, in the last line we take into account the following relations: 7 *S^(— q) r y° = — *S(q) 
and fj(—q) = vKq)l°- On the right-hand side of Eqs. (j3.1j) . f)3.2j) in the expansion of 
induced currents and sources we keep the terms up to the third order in interacting fields 
containing a relevant information on the two-to-two scattering processes. The expansion 
terms of induced currents 3 j A and on the right-hand side of Eq. (j3.1j) in the hard thermal 
loop (HTL) approximation have the following structure: 

j^ a (A, A)(k) = ± g {T a ) bc j*r^ x (k, -k u -k 2 ) A b " \k l ) A cX (k 2 ) 6{k - k x - k 2 )dk 1 dk 2 , 

j* ( ° l2)a (^)(*0 =^ (3.3) 

j^ a (A^m) = 9 2 J ST^ik, -h; q u -qi^ihm-q^q,) 

x 8(k + q\ — ki — q 2 )dkidqidq 2 , 

where (T a ) bc = —if abc . *T flI/ \ is HTL-resummed three-gluon vertex and *r^ (3 is HTL- 
resummed vertex between quark pair and gluon. They represent a sum of bare vertex and 
corresponding HTL-correction. ^^J^^ is an effective (HTL-induced) vertex between 
quark pair and two gluons. The current j A( - 3 ' defines a self-action of soft bosonic field 
considered in Ref. PjJ QB] • The explicit form of HTL-amplitudes in integrands in Eq. ()3.3|) 
can be found in Ref. The superscript (G) points to the fact that in the coordinate 
representation the time argument of external gluon leg incoming in the vertex functions 
and ST™*** is largest. 

Furthermore, the induced sources on the right-hand side of Dirac equation ()3.2|) have 
the following structure: 

Va'^iA V)(g) = 9 (tTj Tjjyfci; ftl -^(fcO^fa) 6(q - gi - fci)d*idgi, (3.4) 

3 Here, we somewhat redefine HTL-induced currents j A< ^ 2 \ j A<,3 \ j*(°' 2 ) and source 7/ 1 ' 1 ) including 
bare vertices in their definition (see, below). 
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V ^(A, A, VO(g) = ^ 9 2 J ST^ih, k 2 - qi , -q^ih^ih)^) 

x 8(q — q\ — ki — k<i)dq\dk\dk 2 . 



Here, *T^ af3 and STj^-^^ are HTL-resummed vertices between quark pair and gluon, and 
quark pair and two gluons, respectively. The superscript (Q) denotes that these vertices in 
the coordinate representation have the largest time argument for external quark incoming 
leg. 

For conjugate induced sources fj^' 1 ^ 0, and fj^' 1 ' 0, by virtue of properties (A. 3) and (A. 4), 
we have from (j3.4j) 

^\A*M-q) = g (t a rj "T«k(-*i; -Qi,Q)^ afl (h)4(-qi) 6(q - ft - fc^M?!, 

(3.5) 

V^\A\ A\ $){-q) = -^ 9 2 l SV^ ji (-q, k u k^(h)A^(k 2 )^(- qi ) 

x 8(q — qi — k\ — k 2 )dqidkidk 2 . 

The 5r^^* ,? and <5r^^ i "' vertices don't exist at tree level and to leading order they 
arise entirely from the HTL In the subsequent discussion we are needed their 

color structure, therefore an explicit form of the vertex functions is given in Appendix A. 

Lastly *V^ u (k) and *S(q) are medium modified gluon (in the covariant gauge) and 
quark propagators 

*V^(k) = -P^k) *A t (k) - Q^{k) *A\k) + £D^(k)A (k), 

*S(q)=h + (q)*A + (q)+h4q)*A4q). 
Here in the first line, the Lorentz matrices are defined by 

Pfiu(k) = g^, v — D^ v (k) — Q fiu (k), Qfj,u(k) = ~— ' Dfiv(k) = ^ 2 , 

A°(k) = — , u^k) = k 2 u^-k^k-u). 

The scalar transverse and longitudinal gluon propagators are *A t,l (k) = l/(k 2 — <5n*''(fc)), 
where SU t,l (k) are transverse and longitudinal soft-gluon self-energies, £ is a gauge fixing 
parameter. Let us assume that we are in a rest frame of a heat bath, so that = 
(1, 0, 0, 0). In the soft-quark propagator the matrix functions h±(q) = (7 =p q • r j)/2 with 
q = q/|q| are the spinor projectors onto eigenstates of helicity and 

*A±(g) = - -= n r ^ . .. , (3.6) 

where 5E±(g) are scalar quark self-energies for normal (+) and plasmino (— ) modes. 
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The system of nonlinear integral equations 1)3.1)1 . (|3.2|) can be perturbatively solved 
(at least in the weak-field limit) by the approximation scheme method. Discarding the 
nonlinear terms in A, if) and ip on the right-hand side of Eqs. 1)3. 1)3.2)1 . we obtain in 
zero approximation 

•^(AOA-CAO = 0, •StfbWfa) = 0, ^(-g) = 0. 

The solutions of these equations denoted by A^ a (k), ip^ l (q) and i>a^(~Q) are the solu- 
tions for free fields. By iterating Eqs. (j3.1|) . ()3.2)) in general case we obtain the series 

oo oo oo 

4i(A0 = £4r 1)a (*). V'i(ff) = E*- 1)< (?) J &(-<?) = E*- 1 ^-?), ( 3 - 7 ) 

s=l s=l s=l 

where A^~^ a (k), ^4* *(?) an d Vq S g) are contributions of g^ 1 order to the soft- 
gluon and soft-quark interacting fields 4 . Substituting expansions ()3.7)l into system 1)3.1)1 . 
()3.2)1 . we obtain iterative solutions of higher order in g: 

to the first order 

A^ a (k) = -*P^(&)[j A ( 2 ^(A^ (3.8) 
to i/ie second order 

Af a {k) = -*P^(A0[j A(3 ^(^^ (3.9) 

and so on. Here, on the right-hand side new effective 5 currents and sources being function- 
al of free fields appear. The effective currents and sources to first order approximation 
()3.8)1 are associated with initial currents ()3.3)1 and sources 13.4)1 . ()3.5)1 by a simple way 

~-A(2)av _ jA(2)au fj\l0) j^P)\ j*(0,2)ai/ /^(0) ^ ^(0) \ = -*(0,2)ai//^(0) ^(0)\ 

4 M)i (A(°)^ (0) ) =^ 1,1)i (i4(°>^ ( ° ) ), ^ 1,1)l '(A*(°)^(°)) S^*^,^). 

4 Note that by virtue of choice of the right-hand sides of Eqs. i(3.1|) and 1(3.2(1 . expansions ((3.7|) accurate 
up to s = 3. For higher s it is necessary to add the subsequent terms of the expansions of induced 
currents j A , and sources r), fj to the right-hand sides of Eqs. 1)3.1(1 . 1(3. 2(1 . 

5 Hcreafter all effective quantities such as currents, sources, vertices etc. will be denoted by the same 
letters as initial quantities with tilde above only. 
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The explicit form of the effective current j A ^ au to second approximation order (J3.9)) was 
obtained in Ref. The remaining effective currents and sources have the following 
structures: 

]* (i, 2 )^(0)^(0) ; ^( 0)) = g ^rf)^ (K - kl - Ql , -q 2 )A^ (h)^(- qi )^\ q2 ) 

x 5(k + qi — ki — q 2 )dkidqidq 2 , (3.10) 
fc^A® A®, ^) = £ j*Tfi:S ij {K k 2 - q u -g)A(°)-^(^)^ ( ° )a2 ^(A; 2 )^(g 1 ) 

x S(q — q± — ki — k 2 )dqidkidk 2 , (3-H) 

X 5(q + qi - q 2 - q 3 )dqidq 2 dq 3 , (3.12) 
and correspondingly, (Dirac) conjugate effective sources are 

x S(q — qi — ki — k 2 )dqidkidk 2l (3.13) 

x 5(q + qi - q 2 - q 3 )dqidq 2 dq 3 . (3.14) 

The effective amplitudes entering into the integrand of the second iteration ()3.1()j) - (|3.14j) 
(and also entering into higher iterations s > 3) represent highly nontrivial combinations 
of the HTL-amplitudes and equilibrium soft-gluon and soft-quark propagators. These 
effective amplitudes can be defined in a direct way starting from the connection between 
the effective and initial currents. Thus for example the effective amplitude in (J3.10)) is 
defined from the relation 

J*(l )2 )a (A (0)^(0)^(0) )(A , )= ^(1,^^(0)^(0)^(0))^) 

+ {/( 2 )^(-*Pj*(°> 2 )(^ 
+ {i *(0,2)a^(l^ 

However, such a direct approach for determination of the explicit form of the effective am- 
plitudes is very complicated even in the second order approximation and as a consequence, 
ineffective. In the forthcoming section we will consider somewhat different approach to 
the calculation of the effective amplitudes, which allows us to avoid many intermediate 
operations and to a certain extent to automate the calculation procedure. This approach 
is extension of the calculating method for purely gluonic effective amplitudes suggested 
in Ref. £Q to the case of an existence of soft-quark degree of freedom in the system. New 
features here, is appearing effective sources 1)3. 12|) . 1)3-1 4)) . which don't have counterparts 
in the expansion of the induced currents and sources in initial system 1)3. lj) . 1)3. 2|) . 
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4 Calculation of effective amplitudes 



The calculating algorithm of the effective amplitudes is based on the following idea. The 
total induced current j^ a +j^ a , and induced source rf a have two representations: by means 
of interacting and free fields that must be equal each other 

oo oo 



E #« a (i4<°>, . . . , A®) + EE2 Ka)a (^ (0) , - - - , AW. g^^ . ^W^W ), 

oo 

rf a [AM = E^* 1)l (A__li^) = ( 42 ) 

s ~ ^ s 

oo oo oo 

E 1)j (a (0) , • • • , a^, v (0) ) + E E ^' 2i+1)i (^ (0) , ■ ■ ■ , ,^,..,^» ). 

Equation (|4.2|) is to be supplemented by a similar equation for conjugate source f\ l a . The 
explicit form of the functions j^ {s)a (A, ...,A) and ~j^ s)a (A^\ A®) was defined in 
Ref. [T] (Eqs. (3.9) and (7.8)). The original currents j^ n ' 2 ^ a (A, . . . , A,ip,ip) and sources 
^(n, . . . , A, t/)) are expressed as 

j*^(A, . . . , A, ^)(k) = i T^^*, - , ft, -ft) 

(n \ n 

fc - E^i + ft ~~ ft Uftdftll^' " = 0,1,..., (4.3) 
»=i / i=i 



X 



A a ^ (h)... A a ^ (k s )^{ qi ) 8[q-q 1 -J2h)dqxl[dk i , s = 1, 2, 



8=1 / j=l 



Here, the HTL-amplitudes *Yf^^ ij for n = and T^;;^/ i for s = 1 represent 
the sum of bare vertices and corresponding HTL-correction. For n > 0, s > 1 we rename 

The interacting fields on the left-hand side of equations (|4.1|) and ()4.2j) are defined by 
the expansion 



A a »{k) = A {0)afl (k) -*V^(k) 



Eif )a (^>-^ (0) ) 



s=2 
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+ E E i* (n ' 2 ^ (0) , - , aw, ^°),v>w,-,^,^ 

1 7 P <■ .. ' •> .. ' 



71=0 1 = 1 



(4.4) 
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E4 S,1) V ( V..,^ (0) ) 



Ls=l 



E E C 1 ' 2l+1)a ( a(0 \---i a(0) ,ip {0 \ip {0 \---,ip 



(oh 



n=0 1=1 



21+1 



Similar expansion holds for conjugate function ip l Q (—q). The calculation of purely gluonic 
effective currents j^ s)a (A^ \ A®) was considered in Ref. pQ. The remaining effective 
currents and sources have the following structure: 



J*(n,20«( A (0) j ;A (0) ? ^(0) ) ^(0) j . ^(0)^(0))^) 



1 1 



+i / *-p(G)aai...a n ,ii...ijji...j ; 



IJ,Hi...fin,ai...aiP 



n! (/!) 2 * 

xA (0)aiw (y • • • A(°^(k n )^(- qi ) . ■4 { S H {-qi)^f 1 n {q'i) ■ --^H 



i 



x 5(k-J2 k i + Efe - Qi) H dk i YidQjdq'j 

i=l ) i=l j=l 



(4.5) 



i=i 



*fi(Q)oi-o s ,ij/, , \ 

1 m...ix s ,ap {Ki,---,Ks,qi,-q) 



si 

x# 9iw (fci) . • • A(°> a ^(fc s )^ )j '(gi) 5( g -gx-E^l^Il 



\ i=i / i=i 

etc. The substitution of Eqs. (|4.4|) . f)4.5|) into (j4.1|) and (|4.2j) turns the last equations into 
identities. Now we functionally differentiate the right- and left-hand sides of equalities 
(jSj) and (|4~2~j) with respect to free fields A^ a>M , ^^(q) and ^i° )i (-'?) considering Eq. (|4~H|) 
for differentiation on the right-hand side and Eqs. (JUSJ), Q4.4j) . ()4.5j) for differentiation on 
the left-hand side, and set A^ a ^ = ip^ 1 = z/>4 = after all calculations. The desired 
effective amplitudes will appear both on the left-hand side and on the right-hand side. 
However principal effect here, is that the effective amplitudes on the left-hand side will 
have at least one less external soft leg than on the right-hand side. This enables us to 
calculate them in a recurrent way. Below we will give some examples. 

The second functional derivative of the induced current j* a with respect to and 
•0(0) yields 



5^f j {q 2 )5^ 



(0)i, 



9 *f £25*(*; ft, -Q2)S(k + q l - q 2 ) (4.6) 



A(°)=V.(°)=i/.( )=0 



g(t a ] 



ij *T^(k;q 1 ,-q 2 )S(k + q 1 -q 2 ) 
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and the second functional derivative of the induced source rf a with respect to and 
•0(0) gives 



g*r(Q^(k-, qi ,-q)5( q - qi -k) (4.7) 



A(°)=V>(°)=^(°)=0 



- q )S( q - qi - k). 

Eqs. ()4.(jj) and ()4.7|1 show that the effective amplitudes of decay of soft-gluon excitation 
into pair of soft quark- ant iquark excitations (or inverse process of pair annihilation into 
soft gluon) coincide with usual HTL-resummed vertex between quark pair and gluon. 
This decay process is kinematically forbidden and hence the 5-functions don't support on 
the mass-shell of plasma excitations. 

The first nontrivial examples arise in calculation of the next derivative. It defines the 
effective amplitudes for the processes with four soft plasma excitations: the process of 
nonlinear interaction of two soft bosonic and two fermionic excitations and the process of 
four soft fermionic interaction. For the first process two effective amplitudes are possible. 
The first amplitude arises from the variation 



5A^^(k 1 )5^ ) \q 2 )8ij [ a ) 
S 2 jM^(k) 5A a >2{k' 2 ) 



a(°)=^(°)=v>( )=o 



5 2 A a >'i{k' 1 ) 



+ 



5A a >'i(k[)5A^(k' 2 ) SA^^ih) 5^ ):i (q 2 )6^ ) \-q 1 ) 



+ 



SA<^(k[)6^,(q' 2 )5^,(-q[) SAi^^ih) 5^ )] {q 2 ) 8$, 



i 0) *( 



-91 



dk^dq^dq'z 



^OZa^Hzi) 8A^^{k{)5^f\ q2 ) <ty. 



■(0)</ 



+ 



dq[dq 2 



d qi dq 2 



Hi(q^\(- q [) Wf j {q 2 ) iA^(fcO#(-ft) 

Here, on the right-hand side we keep the terms different from zero only. Taking into 
account Eqs. (|3.3|) . (|4.4|) . ()4.6|) and ()4.7jl . from the last expression we find the effective 
amplitude determining effective current (j3.10j) 



T 



(G)aai, ij , 



^{G)aa\, ij 



wuap 9i, -92) - (5r m ^ 



(k,-ki, q i,- q2 ) 



[t a ,t a ^ *T^(k, -k + k u -h) *V^'(k - fci) Tj?L(- fi + q 2 ; q u -q 2 ) 
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- (t*T)« Tj^(*; ft, -k - ft) *S 17 ,(k + ft) T«J y/J (*i; -*i - ft) 
+ Tg^fci; -ft, ft - h) *S w (h - q x ) T^(*; -* + ft, -ft), 

where [ , ] denotes a commutator. The second similar effective amplitude results from 
variation 

5^ a [A,^(q) 



6A(°y^{k 1 )5A( ) a ^(k 2 )S'ipf )j {q 



A(°)=i>(V=ip(°)=0 



The calculations analogous to previous ones lead to the effective amplitude determining 
effective source (jH.llj) 



*f fi^iK h; ft, ~q) = Sr^lfik,, k 2 ; ft, -g) 



+ [t ai ,f 



a 2 W *■ 



(4.9) 



-q + q x ; q, -qx) *V UV (h + k 2 ) T^ ww (kx + k 2 , -kx, -k 2 ) 



- {t ai t a2 f *rf£ en {-ki\ q, ~q + h) *S w (q - kx) *r™ 7 ,p(-k 2 ; ft + k 

- (t a H a T Tg?orr(-^ 9> -V + **) * S ^>(<1 - k 2 ) T«J y/J (-* i; ft + fc lf -ft; 



From the definition of this effective amplitude it follows symmetry property with respect 
to permutation of soft-gluon legs useful in the subsequent discussion 



i Mi/i2,0!/3 fc 2, ft, -ft) 



1 /i2Ml,«/3 1^2, fcl, ft, -ft)- 



(4.10) 



The effective amplitude of four soft-quark interaction follows from the variation 

S s V i[A^}(q) 



54f 3 (ft)^if 2 (-ft)^i 1 )n (ft 



J 4(0)=^(0) = ^(0) =0 

5 2 A<^{k[) 



5A^{k'x)5^ l \{q'x) 54f 3 (ft) W n (ft)5V. 



,(o)u 



S 2 V^(q) 



-ft 



dk l dq 1 



^i? 3 (ft)^if 2 (-ft) ^(ft) 



dk! x dq' x 



A(°)=V(°)=V (0) =0 



We kept again the terms different from zero only on the right-hand side. Taking into 
account Eqs. (jH.4j) . (|4.5|) and ()4.7|) . we easily derive the explicit form of the effective 
amplitude in integrand of effective source ()3.12|) 



= (t a r>(t a Y^ T«f(g _ g2; ?2) _ g ) _ q2 ) Tg'(- ft + g 3 ; g l5 -g 3 ) 

+ (tT^T* 2 T«>(? - ft; ft, -g) *2V(9 - ft) *rgf 2 '(-ft + ft; ft, -ft). 
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Here, the following property of antisymmetry with respect to permutation of two last 
soft-quark legs holds 

There is no such a property in general case for permutation of two first soft-quark legs by 
virtue of inequality of the HTL-amplitudes with different time ordering of external legs: 
*F^(q -q x ;q x , -q) ± *T%^(q - q x ; q x , -q). 

For convenience of subsequent references in Appendix B we give the explicit form of 
the effective amplitudes for conjugate effective current j* <I '( 1 > 2 ) a ( y 4*( ), ip^) and Dirac 
conjugate effective sources ^' 1}i (A<°\ A<°\ $(°)) and f} ( ^ ] ^ {0) ). 



5 Boltzmann equation for elastic scattering of soft- 
quark excitations off soft-gluon excitations 

The initial equation for deriving the Boltzmann equation describing a change of the num- 
ber densities of the soft fermion excitations owing to their scattering off the soft boson 
excitations is equation (Eq. (4.4) in Ref. jH]) 

= igjdq'dq x dk x {{^ a {-q)A^{k x )[*Y^ -{k x ;q x ,-q')^{q x )t) 5{q' - q x - k x ) 

- m-q x ) T») fl (-* i; -q x , q)&A*^ 5{q -q x - k x )} 
+ ig 2 j rfg'rfg 1 dA; 1 rfA ;2 {(^(A;0^(fc 2 )C(-g)[5r( < 3) a ''(-A ;i , -k 2 ; q\ -q x )^{q x )t) 

x 5{q' -q x -k x - k 2 ) 
- {A^^kM^^m-q^Yf^^, k x ; -q, q x )t^ a (q')) S(q - q x - k x - k 2 )}. 

Here, ^ = j^q^, the Dirac trace is denoted by Sp, and 5S H (g) is "Hermitian" part of 
the quark self-energy. The angle brackets denote an average with respect to any density 
matrix. 

Using above-obtained expressions (jH.Hj) and for gauge potential and quark wave 
functions by a simple search we extract all sixth-order correlators of free fields responsible 
for the scattering processes of soft fermion excitations off the soft boson excitations to the 
lowest order in the coupling constant. As a guide rule for choice of the relevant terms a 
simple fact will be used. When we make decoupling of the sixth-order correlators in the 
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terms of the pairs (in order to define the product n^N^N^) we keep the terms, which 
correctly reproduce relevant momentum-energy conservation laws. These conservation 
laws encoded in (^-functions of the first two measures of integration in Eq. ([2.7)1 . Besides the 
number of appearing terms can be considerably cut if we note that it is necessary to keep 
only such terms in intermediate expressions, which contain the quark propagators *S(—q) 
and *S(q'). These propagators give later on the terms proportional to Z±(q)5(q° — cu^ ), 
i.e., the factors take into account the existence of soft-quark excitations with wave vector 
q and energies in spite of the fact that the number densities of soft-quark excitations 
riq are explicitly absent on the right-hand side. 

Let us consider the first term on the right-hand side of Eq. (|5.1j) containing the third- 
order correlation function of interacting fields 

ig J dq'dq 1 dk l (^« (-g)A^(*i) *T^ jj ' (h; qi, -g')^'(ft)) " ft " *0- ( 5 - 2 ) 

For this term there exist two substitutions resulting in the sixth-order correlation function 
of free fields with conservation laws (|2.7jl . The first substitution is 

$,(-?) ^vf 1) \A< \A<°\4,M)(-q)% a (-q), 4( qi )^^ )j '( qi ), 
A^ih) — > -*V^(k 1 )j^ a (A ( -°\ A ( -°^)(k 1 ) 

and the second one is 

^(-q)^^' 1) \A<°\A<°\^)*S ia (-q), A^i^^A^ih), 

Here, the product of the correlation functions of free gauge and quark fields appears in 
the following form: 

^4*(0)oia*i ( fcl )^*(0)a 2 M2 ( fca ) A (o)a 3/i3 (Jfe 4 )) (^{-qtyf* (qi)) 

~ { (5°i <I */w«(^) ( 5(^ - fc 4 ) 5 a2a3 i™(A; 2 )5(£; 2 - jfe 3 ) (5.3) 

+ 5 a ^I^ 3 (k l )5(k 1 - k 3 )6 a ^I^(k 2 )6(k 2 - h)}5 i n f3a (q l )6(q l - q), 

where on the right-hand side we use definitions of the correlation functions I^ u {k) and 
T a p(q) in the equilibrium Simple algebraical transformations enables us to lead term 
(|5.2jl to the following form: 

ig 4 f ~(Q)aia2,il 

— J dqidhdk 2 *r mM2 , a/ 3 (h, h; qi, -q) *S aa f[-q) 
x {[t a \ t a *Y l *T^,(q - Ql ; qx, -q) *V W '(h + k 2 ) T^fa + k 2 , -k u -k 2 ) (5.4) 
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- (t a H a *y l Tg^ 7 (-* i; q, -q + h) *S lY (q - h) Tjg y/J ,(-* 2 ; gi + fc 2 , 

- (t a n a ^y l Tg^C-fe; g, -? + fe) ^(g - fc 2 ) T^C-A*; gi + fci, - gi )} 

x T l3 , f3 (q 1 )I^(k 1 )I^(k2)5(q -qi-h- k 2 ). 
In deriving Eq. ()5.4j) we have used property (B.4) for the effective amplitude a p . 

Now we turn to the third term on the right-hand side of Eq. (|5.1jl containing the 
fourth-order correlation function of interacting fields 

ig 2 J dq'dq l dk 1 dk 2 (A*(h)A»(k i )$i(-q)8r^ ii ' (-h, -k 2 ; q', -gi)^'(gi)) (5.5) 

x Kq' - qi-h- fa). 

Here, there exists a unique replacement, which leads to relevant conservation laws ()2.7|) 

-C 1)l (^ ( °U* (0) ^ (0) )^(-?), ^'(gi) -< )5 "(<7i), 

A^(Jfei) A(°^(ifei), A 6 ^(A; 2 ) -> A (0)bi/ (fc 2 ). 
This replacement with decomposition (|5.Hjl results equation (|5.5jl in the following form: 

/ dq x dk x dk 2 *f (fci, fc 2 ; gi, -g) ^T^^t^i, fo; gi, -g) 

x T w (g 1 )J^^(A ;i )J^^(A ;2 )5(g - Ql - h - k 2 ). 

If we now add the last expression to Eq. (|5.4jl and take into account the definition of 
effective amplitude ()4.9|) . then we get 

^- / dqidhdk 2 *f ^"'2 fc 2; gi, -g) *S aQ '(-g) *f J?£5#fai> fc 2; gi, -g) (5.6) 

x T w (g 1 )/ w ^(^i)^ 2 ^(A; 2 )5(g - 9l - fa - fa). 

The remaining second and fourth terms on the right-hand side of Eq. ()5.1|) after similar 
reasoning give in a sum expression complex conjugate to that of ([5.6)1 . To consider the 
simplest interaction process of plasminos with plasmons, it is necessary to fulfill the 
following replacements: 

I^(fa) -> Q^(fa)I l (fa), I^(fa) -> Q^Hfa)l\k 2 ), (5.7) 

T^(gi) - (/i_(q 1 )) w T-(q 1 )5(g° - w" ) + (M^W^-q^O?? + w" ), 

Xa'(-g)^-(^-(q))aa'(*A_(g))*. 

We have determined the function T^/^gi) in the quasiparticle approximation. The spectral 
density T _ (qi) describes the plasmino branch of the soft quark excitations and Y~(— qi) 



19 



describes antiplasmino branch. To take into account weakly inhomogeneous and non- 
stationary in the medium it is sufficient within an accepted accuracy to replace equilibrium 
spectral densities I l (ki), T~(qi) and T~(— qi) by off-equilibrium ones in the Wigner form 
(Eq. (3.6) inRef. 0): 

/'(A*)- 4 = /'(*(, x), * = 1,2, 
T-(qi) — T~ = T-(q 1; x), T-(-q x ) — Tl qi = T^-q^s) 
slowly depending on x — (t, x). 

After replacements ()5.7|) we perform a trivial integration in dq®. For antiplasmino part 
of the function T^(qi) we make a replacement of variable in integrand (|5.6jl : qi — — qi 
!. Finally instead of ()5.6)1 for this special case of nonlinear interaction of soft 



modes in QGP we obtain 

f { 1 

- (i*A_(g))* / dq l dk 1 dk 2 

11 J U z {Ki)1 



9* 



u 2 (k 2 



x 



fi»(* 1 )tt«(fe)*f^;V. k; gi, -g)(Mq)U' 

x^i)^ 2 ) T <g£3$*i, fc 2 ; q 1} -q){h-{(b))w 

x T-^ii^g-ft-fe-fe) (5.* 

x ^(^1)^(^2) Tgr^^x, A; 2 ; -g)(/L(fc))w 



x 



The expression obtained enables us to define the probability of the plasmino-plasmon 
scattering w^~!}Zg(q, qi, fci, ^2) and the probability of pair annihilation into two plasmons 
Wqq^gg I (g, qi] hi, k 2 )- At first we express projectors h-(q) and /i_(qi) in terms of simul- 
taneous eigenspinors of chirality and helicity 

(h-(q))aa' = X] v a(& A)u Q /(q, A), (/i_(qi)) w = J] w^(qi, Ai)^(qi, Ai). (5.9) 

A=± Ai=± 

Furthermore, we take also the boson spectral functions IL and Ii in the form of the 
quasiparticle approximation 

4 = 4M°-4j+'V(^ + <), /^s/'Cki,*), i = l,2. (5.10) 

For the first term in integrand of ()5.8|) in product /L/j. a it is necessary to keep only 
'crossed' terms 
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Employing properties ()4.1Uj) and (B.4) it can be shown that the second term here with 
replacement k\ # k 2 gives a contribution, which is equal to the first term. 

For the second term in integrand of (|5.8J) in product ijLlfc 2 we are needed only to keep 
one term of a 'direct' product 

The 5- functions give us a possibility to perform integration in dk®dk 2 . Now we turn to 
(anti)plasmino and plasmon number densities accordingly setting 

n q = (2 7 r) 3 2Z: 1 (q)T q , l-n- = (2vr) 3 2Z = 1 (q)T q , iV£ = - (2tt) 3 2^Z r 1 (k)4 (5.11) 

where Z_(q) and Z;(p) are residues of the HTL-resummed quark and gluon propagators 
at plasmino and plasmon poles, respectively. 

As was stated above, a presence of the second and fourth terms gives us an expression 
complex conjugate to that of (|5.8|) . It is not difficult to see that in practice this comes to 
a simple replacement of the factor (z*A_(g))* in Eq. ()5.8j) by 

(i *A_(g))* + i *A4q) ~ - 27rZ_(q)%° - w"). 



Taking into account the above-mentioned reasoning, we introduce the following matrix 
element of plasmino-plasmon elastic scattering 

The probability of plasmino-plasmon scattering is defined as 

^^°(q,k i;qi ,k 2 ) = J2 ^'^(-k^k^q^-^^'^C-k^k^q^-q))*. 

AAi = ± 

(5.13) 

Now we consider the problem of deriving the probability for elastic scattering of normal 
quark excitations off soft transverse gluon excitations. In this case instead of replacements 
(|5.7|) it should be used the following those in Eq. ()5.6|) 

/"^(jfci) -> p^ihjPih), i^{k 2 ) -> p^(k 2 )i\k 2 ), 

T W (fc) - (^(qi)) W T + ( qi )5(g° - <) + (^(qi))/3'/ 3 T + (-q 1 )5(g° + <), 
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^ - (Mq)W(*A+(<?)r, 

and instead of (|5.9|) it should be used the expansions 

(Mq)W = 1] Ma(q, A)tt a /(q, A), (/i + (qi))^= ^ ^(qi, Ai)w /3 (qi, Ai). (5.14) 

A=± Ai=± 

Furthermore the transverse projectors (ki) , i — 1, 2 will be also presented in the form 
of expansion in eigenvectors of transverse polarization (Kalashnikov, Klimov [TU] ) 



(5.15) 



= e te)/ % 2 - Af (e^ ■ A*), (A* ■ e(^,&)) =0, z = 1,2, 
where e^*' w are some four- vectors such that e fli (ki, £j) are linearly independent. 

The reasoning similar to that used in deriving the probability of plasmino-plasmon 
scattering (J5.1H)) results in the following expression for desired probability 

*N££>(q, ki; qi, k 2 ) = £ E Tg£ft(-k!, k 2 ; qi , -q)(T~ (-k l9 k 2 ; q x , -q))* 

6,6=1,2 A,A 1= ± 

(5.16) 

where the matrix element for elastic scattering of soft normal quark by soft transverse 
gluon is defined now by the following expression: 



^(-.^...-^^^y^^ 1 ^^) 1 ' 2 ^)" (5 , 7 ) 



— (Q)iii2, ii\ 

%(qi, Ai)T WAt2 Q/3 (-Aij, k 2 ; q u -q)u a (q, A) 



on— shell 



Here, Z + (q) and Z t (p) are residues of the HTL-resummed quark and gluon propagators 
at soft normal quark and soft transverse gluon poles, respectively. 

The calculation of probabilities for other types of elastic scattering processes vJ^^qg^ 
and y^h^g is quite obvious. In both cases it is reduced to corresponding replacements 
of the gluon wave functions in Eqs. 1)5.12)1 . (|5.1Hj) and 1)5.16)1 . (15.17)1 

/Z,(ki)\ 1/2 / u»(ki) \ (Z t {K)\ 1/2 ( e^,&) 



2ul 



i — 1,2 and proper choice of mass-shell conditions on the right-hand side of Eqs. (15.12)1 
and 1)5.17)1 . The factors f (kj) provide renormalization of transverse and longitudinal 
gluon wave functions by thermal effects. 
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The probabilities for more complicated scattering processes, which result in change 
of a type of one or both quasiparticles in a final state (for example, the probabilities 
w qg?Xqg j w^~jl^p and so on) can be obtained by similar replacements of gluon and quark 
wave functions. However, according to remark at the end of section 2 in this case it 
should be performed further research on the existence of solutions for the system of 
energy-momentum conservation laws. 

Matrix elements (j5.12j) and ()5.17|) owing to the structure of effective amplitude (B.3) 
have a simple diagrammatic representation drawn in Fig.^ The first graph in Fig.^ 




Figure 1: The matrix element for soft-quark- soft-gluon scattering. The straight and wave 
lines denote soft-quark and soft-gluon excitations, correspondingly. The blob stands for HTL 
resummation. 

represents a direct interaction of soft-quark modes with soft-gluon modes induced by 
HTL-amplitude ST^ (Eq. (A. 2)). The remaining terms are connected with soft-quark - 
soft-gluon interaction induced by two-quark-gluon HTL-resummed vertex and three-gluon 
HTL-resummed vertex with intermediate virtual oscillations representing s- and t- channel 
contributions. Matrix elements (|5.12|) . ()5.17|) also contain u- channel contribution, which 
it is not drawn in Fig.^ 

Let us consider in more detail a structure of matrix element (|5.12jl . For this purpose 
in the initial definition of effective amplitude (B.3) we proceed from the product of color 
matrices t ai t a2 and t a2 t ai to their symmetric (S) and anti-symmetric (A) combinations 

t a H a2 -> -{t a \t a2 } + -[t a \t a2 ], t a H ai -> -{t ai ,t a2 } - -[t ai ,t a % (5.18) 
2 2 2 2 

where { , } denotes anticommutator in a color space. 

The first term only in the effective amplitude (B.3) is needed some comment. Making 
use of an explicit definition of HTL- induced vertex function (A. 2), we rewrite the vertex 
in the following form (we suppress spinor indices and color indices in the fundamental 



23 



representation): 

Sr^lrKh, fa; ft, ~q) = I {t a \^}5r^f(k u fa; ft, -q)+ l - (fa, fa; ft, -q), 

(5.19) 

where 



qi + ie)(v ■ q + 



■ (<?i + fa) + ie ^ v ■ (qi + k 2 ) + it, 
The functions ST^f^ possess evident properties 

sr ( SS A) (^ **; fc> -q) = ± sr l fi£ A) fa, fa; qi , -g), (5.20) 

W%£ A \-ki, -fa; -qu q) = -f(SrSS A \h, fa; ft, -g))^ . 

Taking into account color decomposition (|5.19j) and change (|5.18jl . we can present the 
matrix element T^ 2 ' Ml as follows: 

T^»' <fa (-k 1 ,k a ;q 1 ,-q) = \ {f>, t"'}"'^ 1 T£>(-k,, k 2 ; q b -q) (5.21) 

+ i[i'«,i<"]« i C^ 1 T(-J!(-k 1 ,k 2 ;q 1 ,-q), 

c„ s {I Cf (c fT ^)}' 



where 

a/2 



x 



tZ w (fa) \ / u^(k 2 ) 



-shell' 



7j(q l5 A 1 )T A ff A f(-k 1 , k 2 ; qi , -q)«(q, A) 
r M (fj 2 (-k 1; k 2 ; q 1; -q) = 6r%lf(k u -fa; -ft, q) (5.23) 



and 



- T«>(fa; -ft - fa, gi ) *S'(-ft - fa) T«)(-fc i; -g, q + fa 

- TgH; -ft + A*, ft) *£(-ft + fa) T»>(fa; -g, g - k 2 ), 

Tfrl (-k lf k 2 ; q 1; -q) = 5T«^(*i. -k; -ft, g) 
+ 2 T«)(g - ft; -g, ft) *P""'(fa - k 2 ) TV^fa - k 2 , -fa, k 2 

- Tjg>(fa; -ft - fa, ft) *S(-ft - k 2 ) T«»(-fa; -g, g + fa) 
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+ T«) (-**; -q x + k u qi ) *S(- gi + h) Tg)(fc 2 ; -q, q - k 2 



The functions possess the properties following from their definitions 

Tg^C-kx.kajqi.-q) = ±Tg^(k 2 , -k l5 q l5 -q). (5.24) 
The advantage of the choice of new color basis ()5.18|) is defined by the following properties: 

({t a \t a2 }{t a2 ,t ai }) ji = 2C F (c F - ^)s ji , 



2N C 

([t a \t a2 ][t a2 ,t a ^ = 2C F (C F + 2^)^'*, (5-25) 

({t ai ,t a2 }[t a2 ,t ai }) ji = 0, 

i.e., crossed terms vanish. Making use of the color algebra (|5.25|h we obtain instead of 
Eq.(jEH 

w qg^qg l) (^ k i> <Ji> k s) = w qfl qg (^ k i 5 Qi > k s) + (q, k x ; q x , k 2 ), (5.26) 

where 



<^(q, k i; qi , k 2 ) = £ | Tg x k 2! qi, -q 

A,A 1= ± 



<«^'-k,.k,:ai.-cri a 

on— shell 



Completely similar decomposition is true and for probability wf^q*' (Eq. (|5.16|l ) and also 
for probabilities wj^ligp, and w^~^~g\ 

The second term in integrand of Eq. ()5.8j) defines the probability for plasmino-antiplas- 
mino annihilation into two plasmons vj^~^gg\ The probability is obtained from ()5.26|) by 
replacement T^f\—ki, k 2 ; q±, —q) — > T^"^(A;i, k 2 ; —q±, —q) (without change of a sign of 
vector qi in spinor t>(qi, Ai)). Accordingly the probability for annihilation of soft normal 
quark-antiquark pair into two soft transverse gluons w^Xgg is obtained by replacement 

T S5ui( -k i> k 2;qi,-q) ~* T ^AAi( k i> k 2;-qi ; -q) also without change of a sign of 
vector qx in spinor w(qi, Ai). 

From Eq. ()5.26|) we see that the processes of elastic scattering and annihilation proceed 
through two physical independent channels determined by a parity of state of two gluon 
system with respect to permutation of external soft-gluon legs (Eq. ()5.24j) ). The symbols 
(S, A) belong to states of two gluons being in even and odd states correspondingly (in 
the c.m.s. of gluons). 

We now return to initial equation (j5.1|) . Let us transform the left-hand side of this 
equation similar to Ref. (Bj. We result in 
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Taking into account above-mentioned, finally we write out the expression for kinetic equa- 
tion defining the change in the colorless soft quark number densities caused by 'sponta- 
neous' processes of soft-quark - soft-gluon elastic scattering and annihilation of soft quark- 
antiquark pair 

V Wfc UX / /i=± b h b 2 =t,l I J 

(5.27) 

+J*i2%?* ) ^ l62) (q, qi; ki, k 2 ) (1 - 

This equation follows from (|2.2j) . ()2.6|) in the limit of a small intensity of soft-quark mode 
ni^ — > and making use of the fact that the boson occupation numbers N^. are much 
more than one, 1 + ~ i = 1, 2. 



6 Probabilities of soft-quark — soft- (ant i) quark elastic 
scattering 

Now we consider the problem of determination of the probabilities for the scattering 
processes of soft-quark excitations off soft-quark and soft-antiquark excitations. For this 
purpose we return to term (|5.2jl and examine the following replacement: 

^^ ' 3) X^ (0) ^ (0) ^ (0) )(-?)^a(-g), ^4 0)i 'fe)> 

A aiAl1 (fci) -> -*^ 1 ^Kfci)j* i (0 ' 2)ai (^ (0) ,V' (0) )(fci) > 

where the effective source r/^ ' 3 ^ is defined by Eq. f|3.14|) . This substitution results in sixth- 
order correlator with respect to free quark fields ip(°> and Furthermore, we decouple 
the sixth-order correlation function in terms of pair correlators (xjj^ijj^) . One can show 
that just two terms of this decoupling reproduce finally necessary factor in integrand: 
$(q + 1i ~ 12 — 93)- Taking into account these terms and the definition of the effective 
vertex *f™ 3 2as , Eq. gJQ , we get instead of Q 

- -^j-Xa'(-g) J dqrfqzdqsSiq + qi - q 2 - q 3 )*r aaia2aa (q, qi , -q 2 , -q 3 ) (6.1) 
x *fiv^ a ' 3 (?,?i, -?2, -g3)T aiQ ' i (gi)T a / a2 (g 2 )T a / a3 (g3). 

At first we will obtain an expression for probability of elastic scattering of plasmino 
off plasmino w(~~^T~)(q, qi; q 2 , q 3 ). For this purpose we set 

*S aa >{-q) - -(Mq)W(*A_(g))*, 
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Ta^gi) -> (^-(qi))^^^^ -w qi ) + (/i + (qi)) aiai T_ qi 5(9 1 +cu qi ) 
T^fe) -> (ft-(4i))a}a i T^<y(g?-w-) + (/i + (4 i )) a { a .Tl^<y(g?+a;-), i = 2, 3. 
Let us especially note that the order of convolution of Dirac indices for soft-quark spec- 
tral density T QlQ /(gi) with the effective amplitudes *T aaia2a3 and *r^,"?^ a , in Eq. (jbUj) is 
the reverse of convolution of Dirac indices of soft-quark spectral densities T a ^ a2 (g 2 ) and 
^0303(^3) (the order of convolution of the last functions coincides with that of convolu- 
tion of the function Y aiQ / (gx) with the effective amplitudes *f| ti ^ 2 Q , /3 and *f^? 1 ^ 2 ,^ in 
Eq. (|5.6jl ). To account for this nontrivial fact in the expansion of the function T aia ' (gi) in 
terms of the spinor projectors in comparison with similar expansion (|5.7j) we introduce new 

functions T qi and T_ qi (with hat above). In the product T aia r (g^T^ « 2 (92) ^0^03 (93) we 
keep only the terms preserved a form of (^-function argument in integrand of Eq. (up 
to a trivial permutations of 4- momenta q±, q2 and g 3 ). There exist three such terms: 

(^-(qi))a 1 a' 1 t qi %? - )(/l^q 2 ))^ Q2 T q2 %£ - W^)(M^))«£«s T i%3 - ^ q3 ) ( 6 - 2 ) 

+ (^+(qi))a lQl T_ qi 5(g° + ^ qi )(/i+(q 2 ))a 2 a 2 T: q2 5(g^ + u^)(h-{q 3 )) a > 3a3 r~ 3 6(ql 
+ (^+(qi)) ai aiT_ qi 5(g? + ^ qi )(^-(q 2 ))a 2 a 2 T q2 5(g^ - cj q2 )(/i + (q 3 ))^ Q 3T: q3 5(g° + 



Let us consider a contribution of the expression in the first line. We substitute this term 
into Eq. (jfi.lj) and perform integration in dqldq^dq®. Further we expand spinor projectors 
/i_(q) in terms of eigenspinors (Eq. ([5.9)1 ) and turn to the (anti)plasmino number densities 
according to rules (j5.11j) supplemented by 

1 - /I" = (2tt) 3 2ZZ 1 (q)f - , n q = (2vr) 3 2Zl 1 (q)f q . (6.3) 

The second term on the right-hand side of Eq. (j5.1|) gives a contribution complex conjugate 
to that of (J5.2)) . The sum of these two terms will differ from (jfi.lj) by a simple replacement 

^(*A_(g))' - \ [(*A_(g))* - A_(g)] ^ 7rZ_(q)%° - u ~). 

Now we introduce the following matrix element of elastic plasmino-plasmino scattering, 
setting 

nTX\(q, qi; -<b, -q») - <? 2 (^f "ft i^r)'* ( 6 - 4 ) 

x « a (q, A)v Q1 (qi, Ai)w a2 (q 2 , A 2 )wa 3 (q3, A 3 )T aaiQ2Q3 (g, g 1; -g 2) -g 3 ). 
The probability of plasmino-plasmino scattering is defined as 

^w^--\q, q i; q 2 , q 3 ) = £ T§SSjq, q l5 -q 2 , -qaXTg^q, qu -q 2 , -q,))*. 

A,Ai,...=± 

(6.5) 
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The diagrammatic interpretation of the different terms of matrix element ()6.4|) is presented 
in Fig.EI 




+ 




Figure 2: Mjzfller-like elastic soft-quark by soft-quark scattering. 



The probability of elastic scattering of plasmino off soft normal quark excitation 

+; 

tions 



Wq q Xq q + ' > is obtained from Eqs. (|6.5|) and (|6.4|) by simple replacement of quark wave func 



I — g — ) ^(Qi^O-^l — - — I u ai (qi,Ai), I — - — I ^ 3 (q3,A 3 )^l — - — I u a ,{q 3 ,\ 3 ), 

and corresponding choice of mass-shell conditions. For deriving the probability of mutu- 
ally elastic scattering of soft normal quark excitations ^qtX'^ it should be also added 
the following replacements to the above-stated those: 



'Z-(q) 



J/2 



v a (q, A) 



'Z+(q) 



J/2 



u a (q, A), 



'Z_(q 2 ) 



vV2 



v a2 (q 2 , A 2 



'Z+(q 2 ) 



vV2 



«a 2 (q2, A 2 ). 



In the remaining part of this section we shall restrict our attention to study of the scat- 
tering processes including only plasmino and antiplasmino. Results for the scattering 
processes with soft normal quark modes follows from previous by above-mentioned re- 
placements of quark wave functions. 

For more detailed analysis of structure of probability ()6.5|) we regroup the terms in 

the effective amplitudes T aQiaaQ , 3 (g, q u -q 2 , -q 3 ) and T^J, ^ (q, q x , -q 2 , -q 3 ). For this 
purpose we replace initial color factors by 'orthogonal' ones 

i i 
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Thus we can represent matrix element (|6.4j) in the following form: 

T^i\(q>qi;-<i2,-q3)= (6.6) 

= I {{tT\t a r 13 + (tT'itT 12 }^ 1 Til| A2A3 (q, q i; -q 2 , -q 3 ) 
+ I {(f )«.(<■)**> - (tT^(tT^}C s l Tg A2A3 (q, q i; -q 2 , -q 3 ), 



where 



4 V N c 

1 aa 1 a 2 a 3 I c I' ( 1i' -q2, -q 3 ) = g <~ StA y — - — j j_X I — 2 — J ( - b - 7 - ) 

x v a (q, A)f Ql (qi, Ai)w Q2 (q 2 , A 2 )t> a3 (q 3 , A 3 ) 
x ' Tg>(- 9 + g 2 ; -q 2 , q) *V vv ,{-q + g 2 ) T^fa - g 3 ; -<7i, <7s) 
± Tgi"(-9 + -gs, q) *V vv ,{-q + g 3 ) T£'( 9l - q 2 ; -q u q 2 )\ 
The functions T AA ' i " A ^ 2A3 possess properties 

T aaV2a 3 (<1> Qi; -Qa, -Qs) = ± T^VaL^ qi; -q 3 , -q 2 ). 
Taking into account a color algebra 

{(t a ) i2i (t a ) i3ii ± (t a ) i2ii (t a ) i3i }{(t b y i2 (t b ) hi3 ± (t b y i3 (t b ) iii2 } 

= 2 {(t a t b y i ti(t a t b ) ± (t a t b t a t b ) ji ) = s ji c F (i t 
[(t a ) i2i (t a f 3ii ± (t a ) i2n (t a ) i3i }{(t & )^ 2 (t 6 ) iii3 =f (t b ) ii3 (t b ) iii2 } = 0, 

we can cast probability (|fi.5j) into a sum of two independent part 

^"^(q, q l5 q 2 , q 3 ) = (q, qu q 2 , q 3 ) + w^l (q, qu qi, q 2 ), (6.* 



"qq^qq 9?— ►ggVMJ M.4-) M^j M«V i ^qq^qq 

where 

«fc£(q, on q 3 ) = £ | Aa (q, qu -q* -q 3 )| 2 on shell - (6.9) 

A,Ai,...=± 

Thus, as in the case of elastic scattering of soft-quark excitations off soft-gluon excitations 
(see previous section), the elastic scattering process of soft-quark excitations off each other 
proceed through two physical independent channels determined by parity of final state of 
soft-quark quasiparticles system. We pay attention to the fact that in the color decom- 
position of matrix element ()6.6|) function T A "^| AaA3 stands with symmetric combination of 
color matrices, and T AA AaAa stands with anti- symmetric combination. 
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The second term in Eq. (|6.2|) defines the process of elastic scattering of plasmino off 
antiplasmino. Reasoning similar to previous one results in expression for the probability 



Wgq^'qq ^(q, qi; q2, qa) analogous to 1)6. B|) . (|6.9J) with replacement 6 of the helical amplitudes 

AAi A2A 3 °y 



'Z-(q) 



1/2 3 

I 5 



'Z_(q< 



1/2 



x 



x v a (q, X)v a2 (q 2 , A 2 )f ai (qi, Ai)w a3 (q 3 , A 3 ) 
rSi"(-g - 915 qi, q) *V vv ,{-q - qi ) Tg£(- ?2 - g 3 ; g 2 , 93) 

± *rg^(-g + g 3 ; -g 3 , q) *V uu ,(-q + g 3 ) Tj$£( gi - g 2 ; <&, -gi^ on _ shdl 

The diagrammatic interpretation of the different terms of matrix element in this case is 
presented in Fig. 01 





Figure 3: Bhabha-like elastic scattering of soft-quark off soft-antiquark excitations. 



Finally, the third term in Eq. ()6.2|) also defines the process of elastic scattering of 

± (S, A) 

plasmino off antiplasmino. The matrix element of this process T^^^^ is associated with 
the preceding matrix element by relation 



(S,A) 



T AA 1 A 2 A 3 (q ! qi; -qa, -qs) = iT^^foqi; -qa, -q2). 

The whole contribution associated with the third term in ()6.2|) can be resulted in previous 
one by a simple replacement q 2 (? 3 . 



^Besides, in the kinetic equation it is necessary to replace (1 — ^- qi ) n q 2 n q3 by (1 — n qi ) 



"qi I '"q2 '"q3 ' 
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Taking into account the aforesaid we write out the final expression for 'spontaneous' 
part of collision terms caused by the processes of soft-quark- soft- (anti) quark elastic scat- 
tering, which should be added to the right-hand side of Eq. (|5.27|) 



E < 

h,f%h=± 



^ ( /i x 9 ; / 2/3) w^ fs \q, qil 03, Os) (l - n^)n^n^ (6.10) 



+ 2jdT q ^ hh) /3) (q,qi;q 2 ,q 3 )(l - n^n<£V°» l 



<\2 q3 



7 Elastic scattering of soft-gluon excitations off soft- 
quark excitations 

Here, we consider derivation of the kinetic equation describing a change of soft gluon 
number densities iV£' 1 caused by their interaction with soft-quark excitations. In this case 
the initial equation is [H] 

9 [ k 2 g^ v - (i + r 1 )*^" + m*» v {k)] dI *} k : x) = (7.i) 



>g_ 

2 



' dk'dk 1 dk 2 {(T b ) cd *T^\k, -k u -k 2 )(A™(k r )A c u (k x )A d x (k 2 )) 5(k - kx - k 2 ) 
+ (T a y d ( *T^ x (k', -k x , -k 2 ))*(Al(k)Al c {k x )Al d (k 2 )) 6{k' -k x - k 2 )} 
+ igj dk'd qi dq 2 {{t b f *T%^(k; q x , -q 2 ) (A; a (k'W a (-q x )^(q 2 )} 5(k + q x - q 2 ) 

- (tT m*V^(k'; gi , -q 2 )h Q U (4(kW a (-q 2 )4( qi )) S(k' + Ql - q 2 )} 
+ ig 2 j dk'dk x dq x dq 2 {^^(k, ~h; q x , -ft) (A^(k')A^(k x W a (-q x )4(q 2 )) 

x 5(k + q x -k x - q 2 ) 

~ (l°[^ )ac (k', -k x ; q x , -q 2 )h°)% (A b ^k)A^(k x )^(-q 2 )^(q x )) 5(k' + q x - k x - q 2 )}. 

In order to obtain the desired equation we must take into account only the substitu- 
tions that lead to the correlation function depending on two pairs of free-quark fields 
^(°)<0(°) ; ip(°)'ip( ) and two free-gluon fields A^°\ A*^ . We consider the first term on the 
right-hand side of Eq. 1)7.1)) containing the third-order correlator of the interacting fields 

- f (T b ) cd J dk'dk x dk 2 *T fluX (k, -k x , -k 2 )(Al a {k')A c v {k x )A{{k 2 )) 8{k -k x - k 2 ). (7.2) 
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The replacement resulting in relevant correlation function has the following form: 

A a/ (ki) - - *V w \k x )j^ 2)a {^\ ^)(h), A d \k 2 ) -> A(°) dA (fc 2 ). 

Here and henceforth, we consider only the replacements that contain gluon propagators 
*Z>«" (A;) and *£> wt (Jfe 7 ). These propagators give later on the factors proportional to 

5(Re *A~ lt ' l (k°, k)) = -Lj Z M (k) [5{u> - ^ l ) + 5(uj + u£% 

The factors take into account an existence of soft transverse gluon modes and plasmons 
with wave vector k and energy (see remark after Eq. (|5.1j) ). 

For this substitution the following product of correlation functions of free fields arises: 

(^-?l)< )i2 (?2)^ 

^ {S tli2 T a2ai (q 2 )6(q 2 - qi)S^T aia3 (q 4 )6(q 4 - q 3 ) 

-5 ill4 T a4ai (g 4 )5(g4 - qi)S l3l2 T a2a3 (q 2 )6(q 2 - g 3 )} ^/""(A*)^! - k 2 ). 

On the right-hand side of the last expression only the second term in braces results 
in collision terms with required conservation laws of energy and momentum, Eq. (|2.9|) . 
Performing the calculations similar to those of sections 5 and 6, we obtain instead of ()7.2j) 

• 4 

_ !|_ [t b^ c] ij * V ^ k )J dkldqidq2 *r^ x (k, —k + h, -h) *V w \k - h) (7.3) 

(Q) — (G)ac, ij 

x * T u>' a >f3>( k - h] qi, -Q2) *IVy,a/3 (*> -h) <h> -92) 
x I X ' x (kx)T aal (q 1 )Tp^(q 2 )S(k + qi-h- q 2 ). 

For term (|7.2jl there exists another relevant substitution 

A* a »(k') -> - (*2? w '(jfcO)*V* (1 ' 2)a (A* (0) ,^ (0) , ^ (0) )(fc'), ^(h) -> A(°) a '(jt 1 ) 1 

^(fe) - -*P AA '(A; 2 )j A r ' 2)d (^( ),^( ))(A ;2 ). 

The expression obtained for this replacement by a simple change of integration variable 
can be resulted in (|7.3|) . 

Furthermore, we consider the third term on the right-hand side of basic equation (j7.1j) 
also containing the third order correlation function of interacting fields 

ig (t b fjdk'd qi dq 2 Tg^*; q x , -q 2 ) (A; a (k'W a (- qi )^(q 2 )) 8(k + q x - q 2 ). (7.4) 
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For this term there exist two substitutions leading to sixth-order correlation function of 
free fields with conservation laws ([2.9)1 : 

^(ft)^-^(fe)^ 1,iy (^^^ (0) )(fe), 

and 

A*"^*/) -> -(*P^'(^))*^* (1 « 2)a (A*W,^°J,V' (0) )(*')» 

^ €' 1)l (^ (0) ^ (0) ) *S a , a (- qi ), ^(q 2 ) -> ^(g 2 ). 

These replacements after simple algebraic transformations lead term (|7.4j) to the following 
expression: 

-ig 4 *V^'{-k)Jdk 1 dq 1 dq 2 (7.5) 
x | ( ^ C )ii Tj^ 7 (*; 9l , -* - 9 i) *S 1 y(k + qi ) T^CAn; g 2) -fc - g 2 ) 
+ W *ri%#i; -gi, gi - fci) *S 11 >{k - g 2 ) T^(*; -* + g 2 , -g 2 )} 

~(G)ac,ij 

x * r M 'A', Q /3 ?i, -0z)J (^i)T Qa '(gi)T^ /3 (g 2 ) <J(fc + g x - fex - g 2 ). 

Finally, we will consider the terms containing fourth-order correlation functions on the 
right-hand side of Eq. (|7.1|) . Here, there exists a unique relevant substitution 

A*^(k') -> - (*X> w '(jfc , ))%** (1 ' 2)a (A* (0) ,^ (0) , V (0) )(^) ; ^(Jfci) -> v4 (0)ci/ (A:i), 

This replacement results the first term containing the fourth-order correlation function of 
interacting fields in the following form: 

- ig\ *V^\k)Yj dhd qi dq 2 6T% )b :;;i(k, -h; q u -ft) ^^(k, ~h; qi, -g 2 ) (7.6) 

x I X ' X {ki)T aa/ (q 1 )Tf 3 ^(q 2 )5{k + q 1 - h - q 2 ). 

Now we add together expressions obtained (|7.Hjl . (|7.5jl and (|7.fj|l . Taking into account the 
definition of effective vertex (|4.8jl . we find that the sum of the first, third and fifth terms 
on the right-hand side of Eq. (j7.1|) can be presented in the form similar to Eq. ()5.6j) 

- tg 4 *V^\-k) J dhd qi dq 2 *f ^rj(fc, -h; q h -g 2 ) *f ( ffi$(k, ~h; g x , -g 2 ) (7.7) 

x J A ' A (A;i)T C(Q /(gi)T ) g/ ( 9(g 2 ) + gi - &i - g 2 ). 
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Further, we transform expression ()7.7|) within the scheme suggested in section 5. Let us 
obtain at first the probability of plasmon-plasmino elastic scattering. For this purpose 
at the beginning we extract a purely plasmon part from gluon spectral density I x ' x (ki) 
setting 

7 V A(fci) _ u x \h)u\h) ^ m _ <} + ^ m + ^ 

( *V^\k)Y -> - ^( fc K ( fc ) ( *A'(fc))*. 

By virtue of 5-functions we can perform integration with respect to dk\. In term containing 
li kl we make a replacement of integration variable: kx — ► — k x (wL — > wL)- Further, we 
extract the plasmino part of the functions T aa i(qi) and Ypip(qz) 

T aa f(q 1 ) -> (/i-(qi)) aQ 't ~ <J(gJ - w~ ) + (/i+(qi)) Q a'T_ qi 5(g 1 + w~ ), 

T w (g 2 ) -> (/i-(q 2 ))/3' ) 8Tq 2 5(g^ - w~ ) + (/i+(q 2 )) w T: q2 5(g° + w~ ). 

In the product Y aa >(qi)Tp/p(q2) for collision term containing conservation law S(k + qi — 
ki ~ Q2) (Eq. (|7.7j) ) it is necessary to keep only two terms of a 'direct' production, i.e., 

T aa /(gi)T w (g 2 ) -> (/i_(qi)) aa /f ~ - cj~ )(/i_(q 2 )) w T q2 - w~ ) (7.8) 

+ (/i + ( qi ))^f I qi 5( g ? + c. q j(Mq 2 ) WT: q2 % 2 ° + ^ q2 ) 

and in the term containing conservation law 5(k + g x + k\ — q 2 ) it is necessary to keep 
only 'crossed' term 

T aa /(gi)T w (g 2 ) -> (/i + (qi)) QQ /f _ qi %? + ^ qi )(Mq2)WT q2 <5(g 2 - w~ ). (7.9) 

The first term on the right-hand side of Eq. (J7.8)) is associated with the scattering process 
of plasmon off plasmino and the second one is associated with the scattering of plasmon 
off antiplasmino. Finally, contribution (|7.9|) defines the process of creation of plasmino- 
antiplasmino pair by two plasmon fusion. 

Let us consider for the sake of definiteness the first term on the right-hand side of 
Eq. (|7.8|) . We substitute this term into Eq. (|7.7|) and perform integration in dq^dq®. For 
spinor projectors /i-(qi) and fo-(q 2 ) we use representation (|5.9|) . The spectral densities 
^ki> ^qi; an d Y~ 2 are replaced by plasmon and (anti)plasmino number densities according 
to rules (j5.11|) and ([6.3)1 . The remaining conjugate terms on the right-hand side of Eq. (17.11) 
lead to the expression similar to ()7.7j) with a unique replacement ( (k))* -> *£> w (fc). 
Thus in the case of plasmon-plasmino interaction for a weak-absorption medium when 
Im*A _1 '(fc) — > 0, this results in simple replacement (*A z (/c))* by 

(*A l (jfe))* - *A'(A;) = -2zIm*A*(A;) ~ 2«sign(A; ) (^pl <J(*° ~ 4)- 
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Here, on the right-hand side we keep only positively-frequency part. 

Taking into account above-mentioned, we introduce the following matrix element of 
plasmon-plasmino scattering: 

M S - (k ,- ki;qi ,- q2) ,^^y7^y / 75#iy / 75#)y /2 ( , 10) 



\H; v 2< / v - / v - 

~ (G)aai, i\%2 

vp(q2, A 2 ) *r wiia/3 (k, -ki; q u -q 2 )v a (qi, K 



on— shell 



The probability of Compton-like plasmon-plasmino scattering in this case is defined as 
(5 a6 «fcy(k,q 1 ;k lj q a ) = £ ^ lfe (k,-k x ;q 1 ,-q a )(]<^(k I -k 1 ;q 1 ,-q a ))*. 

Ai,A 2 =± 

(7.11) 

The different terms of matrix element (|7.1U|) have a similar diagrammatic interpretation as 
is depicted in Fig. ^ with relevant replacements of momenta and HTL-resummed vertices 
only. The probability for the elastic scattering of soft transverse gluon excitation off soft 
normal quark excitation is defined by the following expression: 

S ab w^(K Oil k Xj q 2 ) = E E M !£Tj 2 (k, -ki; qi , -qa)(M*£™ (k, -k l5 q x , -q 2 ))*, 

£,£l=1,2 Ai,A 2 =± 

(7.12) 

where now the matrix element is 
M, a ; AlA2 (k,-k i;qi ,- q2 ) = , (^-| (7.13) 



~ (G)aai, i\%2 

Ufs(q2, A 2 ) T wlja/3 (fc, -A*; g 1? -g 2 )w Q (qi, AiJ 



on— shell 



The probabilities for other types of the elastic scattering v^T^Z', t&^ra are obtained 
from previous by the appropriate replacements of gluon and quark wave functions and 
choice of mass-shell conditions on the right-hand sides of Eqs. (|7.10|) and (|7.13j) . 

As in section 5 we regroup the terms in the effective amplitudes *T a g and 
*f^^g* 1 * 2 according to transition ([5. 18)1 . We present the HTL-induced vertex 5T^^ ai 
in the form similar to (|5.19jl . where now by virtue of initial definition (A.l) instead of 
cT^ 5 '" 4 ) we imply 

W1 J 47r (v ■ qi — it)(v ■ q 2 + it) 

1 

± 



v ■ {qi + k) + it v ■ (qi — ki) — it) 
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Taking into account this expression, we present matrix element (j7.10|) in the following 
form: 



M aa,M 2(k) _ ki . qi; _ q2) = -{t a \t a }^Cs 1 Mf 1 l(k 7 -'k 1 ;C ll ,- q2 ) 

+ ^[t ai ,t a ] i2il 6 1 Mif{ 2 (k,-k 1 ;q 1 ,-q 2 ), 



(7.14) 



where 



M ( 5 ,^ v . s 2 ^ /Z i (k)\ 1/ Yz i (k 1 )Y / Y Z_(q 1 ) y / Y Z-(q2) \ 

MU 2 ; (k,-ki;qi,-q2) = g C s ,a — r hrr 



1/2 



V2, 



V2, 



1/2 



24 ; V 2^ 



ki 



2 / V 2 / 



u"{k) \ ( u^ih) 



(7.15) 



q 2 , X 2 )Mf f ; i A) {k, -ki; qi, -q 2 )^(qi, Ai) 



The functions -M^y^ are defined by the following expressions: 

Mg> (k, -k i; q 1; -q 2 ) = 5Tf^ s \-k, h; -ft, ft) 

- Tg)(-*i; fci + ft) - ft) *Tf\-k; k + Ql , -ft) 
+ *Tf\-k; k - q 2 , q 2 ) *S(-h + ft) T(«K; ft, -ft + h) 



on— shell 



and 



M^\(k, -k x ; qi, -q 2 ) = arjg^C-*, ?2 ) 
+ 2 Tf - ft; -ft, q 2 ) *V W ' \-k + h) T^-fc, k u k- h) 

- Tg>(-*i; -ft, fci + ft) - ft) Tj?>(-*; fc + ft, -ft) 

- Tf >(-*; fc - ft, ft) *S(—ki + ft) T«)(-fc i; ft, -ft + fci). 
Here, Tp is index of the fundamental representation. 

By virtue of representation ()7.14j) probability (|7.11J) can be introduced in the form of 
a sum of two independent parts similar to (|5.2fij) 



it', 



g q ^!gq(K Qi ! ki , q 2 ) = w^ gq (k, qi ; ki, q 2 ) + «^ w (k, q l5 ki, q 2 ), 



91^91 



91^91^ 



(7.16) 



where 

to J g^(k l q 1 ;k ll q 2 )= £ | M £f (k, -k l5 qx, -q 2 ) 

Ai, A 2 =± 

In deriving this expression a color algebra was used 



on— shell 



tr ({t ai , t a }{t\ t^}) = 25 ab T F (C F - JL) . 
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tr([t a \t a }[t b ,t ai }) = 2S ab T F [C F + 



2N r 



tr({t a \t a }[t\t ai }) = 0. 

The decomposition (|7.1fij) is true also for probability (|7.12jl and for other probabilities of 
elastic scattering: i^q^gq and w^jj^. 

Now we consider the second term on the right-hand side of Eq. (|7.8j) . As was men- 
tioned above this term defines the process of plasmon scattering off antiplasmino. The 
reasoning similar to previous one results in scattering probability Wg^^q(k, ki; qi, q2), as 
was defined by Eq. (|7.1fi|) . where now 



(S,A) 
gg^gq 



(k,qi;ki,q 2 ) = J2 \ M A!A 2 ( k > ~ k i; Qi> -Qa] 

Ai,A 2 =± 



on— shell 



and 



M^(k, -k i;qi ,-q 2 )=^4^' 



A1A2 



' Zi(k) 



V2/r7 /l \\l/2/r7 / \\l/2/ v / \ \ 1/2 

\ /Z/(ki)\ /Z_(qi)\ /Z_(q 2 )\ 



2^ J \ 2 / V 2 



(7.17) 



(qi, AOA^gt^Ck, -ki; -q 2 , qiMq 2 , A 2 ) 



on— shell 



In the kinetic equation it is necessary to replace n qi , n q2 by n qi , n q2 . The integration 
measure here remains invariable, as it was defined on the first line of Eq. (|2.9|) . The 
diagrammatic interpretation of the different terms entering into matrix element (|7.17|) is 
presented in Fig. El 





Figure 4: Compton scattering of soft-gluon off soft-antiquark excitations. 

Finally we will consider the contribution of term (|7.9J) determining a pair produc- 
tion gg — > gg. The probability of this process WggL+qq (k, ki; qi, q 2 ) is defined by an 
expression, which is similar to Eqs. (|7.1fi|) . (|7.15|) with a unique difference 7 that functions 

7 In the kinetic equation we replace 1 — by and 1 + iV^ by Af£ . The integration measure is 
defined by the second equation in Ij2.9|) . 
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.MjJ'^k, -ki;q!,-qa) in (TTToT) should be replaced by Mfy^ (k, k x ; -q 1; -q 2 ). The 
diagrammatic interpretation of the different terms entering into the matrix element of 
pair production is presented in Fig. [HI 




Figure 5: Production of soft quark-antiquark pair. 

Now we return to initial equation (|7.1|) . We transform the left-hand side similar to 
[THj . As result we have 

Taking into account above-mentioned, we write out the final expression for the kinetic 
equation defining the change of the colorless soft-gluon number densities caused by 'spon- 
taneous' processes of soft-gluon -soft-quark elastic scattering and pair production by fu- 
sion of two soft-gluon excitations 

(¥ +v ^*T-£ •£ {/^.«^«(k, qi ;k 1 , q2 )<'( 1 -< ) )<' ) 

\ / h=t,l /i,/ 2 =±l 

+jdr^ «#*£ A >(k, q i; k 1; q 2 )<) (l - n£>)n£> 

+ / dl£%™ ^ l/2) (k, k i; q 1; q 2 )< l) n^)ng)}. 

This equation follows from (j2.3jl . ()2.8)1 in the limit of a small intensity of soft-gluon 
excitations iV^ — > and making use of the condition 1 + JVg ~ iv£ l} . 
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8 Generalized formula of emitted radiant power in 
QGP 



In previous sections we have explicitly shown how one can define the probabilities of 
the simplest processes of nonlinear interaction between the soft collective modes Fermi 
and Bose statistics in QGP. This has been done by immediate extracting all possible 
contributions to matrix elements of these processes. However the use of such a direct 
approach in the general case of the processes including an arbitrary number of soft fermion 
and boson excitations becomes ineffective because of awkwardness and complexity of 
computations. In this case the method developed by Tsytovich [201 E] in the theory 
of usual plasma known as the correspondence principle is more convenient for deriving 
an explicit form of the scattering probabilities. We have already used this approach for 
deriving the general expression for the probability of 2n + 2 -colorless soft-gluon decay 
processes, n — 1,2, ... pQ. Here, we don't set ourselves a purpose to obtain similar general 
expressions for the probabilities of the processes including an arbitrary (even) number 
of soft boson and fermion excitations. We concentrate our attention on deriving the 
most general expression for the emitted radiant power, which would take into account an 
existence of fermionic sector of plasma excitations in QGP. Having at hand this expression 
and following standard scheme stated in section 4 of the paper pQ with employing general 
expressions for effective currents and sources (|4.5)l . it is not difficult to obtain the general 
expressions for the probabilities in this more complicated case. Thus the problem stated is 
reduced to calculation of explicit form of the effective amplitudes in integrand of effective 
currents and sources ()4.5|) . Calculating algorithm of the effective amplitudes was proposed 
in section 4 of the present paper. 

For derivation of the probabilities of soft-gluon decay processes in [l] we have used the 
following formula for the emitted radiant power: 



where V is a spatial volume of integration. Averaging over the time is made for elimination 
of oscillating terms in I. Thus the emitted radiant power is defined as a work per time 
done by radiation field with color current (in this case j a (x, t) = j Aa (x, t)) creating this 
field. The chromoelectric field E a (x, t) is associated with j a (x, t) by the field equation. 

We propose that generalization of the above-written formula for emitted radiant power 
taking into account an existence of soft fermion excitations in the system is 




r/2 



// 



dxdt (E a (x,t) -j a (x,t)), 



t/2V 



X — 1b + 1-Ti 
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where 



r/2 r/2 

1 B = lim — / /cM(E a (x,t)-j Aa (x,t)> + lim — / f dxdt (E a (x, t) ■ f a (x, t)) 

t, V^oo TV J J t, V^oo TV J J 

— r/2 V -t/2V 

(8-1) 

= lim £^f- fdkdk°(E a {k,k°)-j Aa {k,k°)) + lim fdkdk (E a Ck,k°)-f a Ck,k )) 

t,V^oo TV J t,V^oo TV J 

and 

r/2 

X^ = -I Ijm^ — / |rfxdt{((^(x,t)^r^(x,t)) + (^(x,t)0 o ^(x,t)r)} (8.2) 

-t/2V 

= -1^^ JdqdqVq^m-q)^)) - (^(-g)^(g))} 

+ f r,feo / dS^f W<^(-9K(*)^(9')> S(q -q'-k) 

~ (n i HVa(-q')A a (k'mq)) 6(q -q' + k)}. 

Here, on the first line 8 of Eq. (|8.2j) we make use of covariant derivative D = d/dt + 
igA^x, t)t a . For Fourier transformation we have used a relation 

r/2 

lim y e iu; *dt = 2ttS(uj). 

-r/2 

The appearance of the second term on the right-hand side of Eq. (18. lj) containing the 
induced current j* a (x, t) is rather evident. The appearance of contribution ()8.2|) is a 
principal new feature, and therefore our further attention will be concentrated on a proof 
of the fact that expression ()8.2|) correctly reproduces the scattering probabilities obtained 
in sections 5 and 6 within the framework of the generalized correspondence principle. 

According to the correspondence principle for determination of the scattering proba- 
bilities of the processes of a type 

qg -»• Q9, qq~> qq, qq -> qq, etc 

8 The authors are grateful to reviewer pointed to the necessity of entering covariant derivative -Do 
instead of usual time derivative d/dt in Eq. H8.2JI that is necessary for maintenance of gauge symmetry. 
The terms in the emitted radiant power Xjr including A§ should be taken into account in the consideration 
of more complicated scattering processes than the processes presented in this work. For this reason in 
the subsequent discission these terms will be omitted (or one can assume that we work in Ao-gauge). For 
recovery of Lorentz symmetry of Eq. Ij8.2|l instead of Dq it should be used (u-D), where u is four- velocity 
of plasma. 
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the expression (|8.2|) should be compared with the expression determining a change of 
energy of soft fermionic plasma excitations caused by the spontaneous processes of soft- 
quark and soft-gluon waves emission only 




In deriving the right-hand side of ()8.3|) we have used equations (I5.27|) and (J6.10|) . To 
proof ()8.2|) it is sufficient restored some scattering probability on the right-hand side of 
Eq. (J8.3)) . For defmiteness we consider the probability of elastic scattering of plasmino off 
plasmon entering into the first term. 

At first step we transform Eq. (j8.2j) to more suitable form. The Fourier components of 
quark fields ip l a (q), ^a(-q) are associated with induced sources rf a (q) and fj l a (—q) by the 
Dirac equation 

V&Gz) = - •Sapto) vfa), &(-<7) = *s Pa (- q ). 

Substituting the last expressions into Eq. ()8.2|) . we obtain instead of (|8.2j) 

X F=- l 2 T $™J^f- /^% Q (t(-?)^(?)){*^(-?)+^(?)}- (M 

By virtue of the properties of the quark scalar propagators *A±(—q) = (*A =F (g))* the sum 
in braces in integrand of ()8.4j) can be written in the form 

2zIm(*A+( g ))(/i + (q))^ + 2zIm(*A_( g ))(/i_(q))^. 

Furthermore, we express the projection operators h±(q) in terms of simultaneous eigen- 
spinors of chirality and helicity 

0+(q))a/3 = M «(q> A )M4 A )- ( /i -(q))a / 3 = v <*(q> A )> 

A=± A=± 
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and approximate imaginary parts of the scalar propagators for a weak-absorption medium 
by the following expressions: 

Im *A ± (g) ~ T vr Z±(q) 6{q° - wj) ± vr Z T (q) + wj). 

It is not difficult to see that the term containing function S(q° + tuT) after substitution 
into (|8.4p . performing integration in <igo an d replacement of variable q — > — q results in 
contribution to emitted radiant power Tjr exactly equal the contribution of term with 
S(q° — Wq). Taking into account above-mentioned, we define the final expression for the 
emitted radiant power Xjr, which is more convenient for concrete applications 

vi 



2> = 2tt lim^ i-^- £ Jdq {< Z + (q) (| t(-<?K(q, A)| J,^) (8.5) 

+ ^ q Z_(q)(|^(-gK(q,A)|; 0=Ljq )}. 



For derivation of the probability for the elastic scattering of plasmino off plasmon in 
the last term on the right-hand side of Eq. ()8.5|) we perform the following replacement: 

where the effective source Va' 1 ^\~ ( l) is defined by equation (jH.lHjl . As a result of such a 
replacement for 'plasmino' part of expression ()8.5j) we have 

(X^)piasmino = *g ^ hrn^ — — J dq cu q Z_(q)(T mM2ja7 (An, k 2 ; q h -q)v a (q, A 



x TOg^Cfci, *£; q[, -?)fy(q, A)) (A*W ai w(^)^* (0)a2M2 (fe)^ (0)a ' i ^(A ; i)A( ) a ^HA ; y) 

x <^ 0)il (-<2l)V>y K (<Zi)> % -qi-h- h) 5(q - q[ - k[ - k' 2 ) dq x dk x dk 2 dq' x dk\dk' 2 . 

Here, the product of correlation functions is transformed by formula (|5.Hjl . Performing 
integration in dq^dk'^dk'^ we obtain instead of the last expression 



(X^plasmino = V hm ^£ f dq Wq Z_ (q) ( ( fc a ) T y7 (ft ) (8.6) 



(2tt 

A=± 

('^^fa. k * 9i, -g)^a(q, A)) (Tg;^^, k 2 ; ft, -g)^(q, A 



x 

2 



x [5(q — gi — /ci — /c 2 )] dq\dk x dk 2 . 
By <5-f unction squared we mean 

[% - q x - k x - k 2 )} 2 = — y rV5(q -q x -k x - k 2 ). 
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In functions I fJ ' 1,J ' 1 (ki), I^ 2 ^' 2 (k 2 ) and Ty 1 (qi) we extract plasmon and plasmino part 9 
according to Eq. ([5.7)1 . Further, going to the plasmino and plasmon number densities 
(Eq. (|5.11jl ). we result equation ()8.6|) in the final form 

,rr \ [ dq dq± dki dk 2 , s a c , n r , / _ / 

(X^) p i asmino = 2 y (27r) 3 (27r)3(27r)3(27r) 3 (2tt) 5(q + kx - qi - k 2 ) <5(cu q + w kl - Wqj - 

x w- f E ^'^(-ki, k 2 ; q 1; -q) (TX' in (-ki, k 2 ; q 1; -q))*) < < , 

\A,Ai=± / 

where the function T"\" 2,ni is defined by Eq. (|5.12|) . The expression obtained exactly 
coincides with corresponding term on the right-hand side of Eq. (|8.3j) . 

Thus suggested expression ()8.2|) correctly reproduces the scattering probabilities ob- 
tained by a different way. Formulae for emitted radiant power ()8.1j) and ()8.2|) open a 
way for more direct and economic deriving the scattering probabilities of the other higher 
processes of nonlinear interaction of soft fermion and boson excitations in QGP. 



9 Linearized version of Boltzmann equation 

In this section we briefly consider the linearization procedure of self-consistent system 
of kinetic equations ()2.2|) and (|2.Hj) . We examine plasmino kinetic equation ()2.2j) with 
generalized rates ()2.5j) and ()2.6|) , where we will consider only the contribution with f\ = — 
and bi = b 2 = I. We assume that the off-equilibrium fluctuations are perturbative small 
and present the number densities of colorless (anti)plasminos and plasmons in the form 

n~ = n; q (q) + Sn~, n q = n~ (q) + 5n^ N\. = N l eq (k) + SN^ 

where n~(q) = n~(q) = (e UJ i/ T + and N l eq (k) = (e UJ k/ T — l) -1 . Parametrizating 
off-equilibrium fluctuations of the occupation numbers Sn~, 5n~ and 5iV k as follows 

x - _ dn- jq) _ 1 _ _ _ 

q= du^ Wcl ~ T e Q^( ~ n eq\ C V) W (l > 

dw-(q) 1 



^ e(? (q) 

_ dn-(q) L 

q = CLJ^ WCI = T Ue *( q )( 1 ~ n eq\ C V) W c l > 



q 



^ fj^K = ^<(k)(l + <(k))W k , 

9 In the product of the spectral densities JjL considering quasiparticle approximation (|5.10(l we keep 
only the terms corresponding to the scattering of plasmino off plasmon. 
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we derive from Eqs. (|2.2|) . ()2.5|) and ()2.6|) . after simple algebraic transformations the lin- 
earized Boltzmann equation for the function w~ 

dw~ dw~ 

dt + q ax 1 j 

o/>r n-( qi )iVi q (k 2 )(l + iVi ^kQ) , _ , 

~ 2 / ^ — — 7-x 2 «W-«(q, ki; qi, k 2 )(w; q + W kl - u> qi - W k2 } 

J '''eg v 4 / 

- / ^™ - — — 73 w gg-- 99 (q> qi; k i> k 2){^ q + - w kl - w k2 | 

U (1 - w- (qi))w- (q 2 )7i- (q 3 ) _ _ _ 

- / dT qq ^ qq — — Wgg->gg(q, qi| ki, k 2 ) {u> q + W qi - W q2 - W q3 } 

J n eq\H) 

of^T (1 - w~ (qi))n~ (q2)n- (q 3 ) _ _ 

-2 dT qq ^ qq a —J- ^ Wgg-OT(q,qi;ki,k 2 ) {w q +w qi -w q2 -%3}- 

J ^egiqj 

In deriving Eq. (|9.1|) we have used the identities 

(1 - n" (q))n- ( qi )(l + A^ki))]^) = n" (q)(l - n~ (q 1 )X(k 1 )(l + <(k 2 )), 

(1 - n~ (q))(l - (qOXlk^iV^k,) = n~ (q)n~ (q x ) (1 + <(k 1 ))(l + <(k 2 )), 
(1 - n~ (q))(l - n~ (qi)) ™~ (q2)n~ (q 3 ) = n" (q)n~ (qi) (1 - n~ (q 2 ))(l - n~ (q 3 )), 
which hold by virtue of corresponding conservation laws of energy. 

In the remaining part of this section we make a little digression from the general line 
of our present work and consider the problem of the minimal generalization of Boltzmann 
equation ()2.2j) to the case of the dynamics of color soft-quark and soft-gluon excitations. 
For this purpose we assume that a time-space dependent external perturbation (e.g. exter- 
nal color current j^ ta {x)) starts acting on the system. In the presence of the external color 
perturbation the soft gauge field develops an expectation value (A®(x)) = A^(x) ^ (but 
nevertheless 10 (i/; l a (x)) = = 0) an d the soft-quark and soft-gluon number densities 

acquire a non-diagonal color structure. In this case, we expect the time-space evolution of 
n q = (n q u ) to be described by the (matrix) plasmino Vlasov-Boltzmann equation instead 
of (JOJ) 

(V t + v q • P x )n q - l -g{ (g(x) + (v- x B(x))) ., V Pi n q } = - C [n~ N l k ], (9.2) 
where T>^ is the covariant derivative acting as 

T>fj,n~ = <9 M n q + ig[Ap(x),n~], A^{x) = t a Al{x). 



10 We can formally suggest that in the system there exists also an external color source r;® xtz (x) produced 
mean quark field (ip l a (x)) ^ (see, for example [22J ) - Here, we don't consider such a case for simplification 
of the problem. 
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The functions S l (x) and B l (x) are mean chromoelectric and chrome-magnetic fields. The 
collision term C [n~, iV£] has the following structure: 

C K,N l p ] = \ {n-,r^[n- JV*]} n~), ifVq, Aft} - . . . , (9.3) 

where ri _) [n q ,iV^ = (if [n~ , iV£]) and r[ -) [n" , iV£] = (rf [n q ,iV^) represent the gen- 
eralized decay and regenerating rates of color plasminos, respectively. The dots on the 
right-hand side of Eq. ()9.3j) is refereed to possible contributions containing commutators 
of [Re£#,n q ] type, where £_r is retarded soft- quark self-energy. 

As above, we restrict our attention to the process of the elastic scattering of color 
plasmino off color plasmon. In this particular case the decay and regenerating rates have 
the following form: 

rf KX] = E fe^°T;^ 

A,Ai = ± J 

x (1 - n-y^iNlJ^il + Nl 2 ) a 'i a \ (9.4) 
rf'K'^k] = E fe^T;^ 

A, Ai = ± 

x (n- ) i i il (l + A^J< 01 (iV£ 2 ) 02a 3, 

where matrix element x^" 2,Ml is defined by Eq. (|5.12jl . We recall that this matrix element 
has a physical sensible color decomposition: 

T---(-ki,k 2 ; qi ,-q) ^r,i" i rC s 1 TS(-k ll k 2 ;q ll -q) (9.5) 

+ i[t-,t-r^ 1 T^(-k 1 ,k 2 ;q 1 ,-q). 

Here, helical amplitudes T ( ^ A) are given by Eq. lET22t and C s ,a={Cf(C f tI/(2N c ))/2} 1 / 2 

Furthermore, we consider the linearized version of Vlasov-Boltzmann equation ()9.2|) . 
We write the number densities of color plasminos and plasmons as follows: 



dn ( o ) 



(iV£) af> = iV^(k)5 ab -^i W^ c (T c ) ab . 

After some cumbersome algebraic transformations we derive the linearized kinetic equa- 
tion for color plasminos from Eqs. ()9.2|) - ()9.5|) 

V t + v q • X? x ) Wq = - g (v- • - 5C [u> q , Vt^], (9.6) 
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where plasmino-plasmon linearized collision term 5C has the form: 

*c K ,XH E "'° (q ' )Af '°' k2) ' 1 . + W [{ixgr + iTinn> q - 

A,Ai=± J n eq\H) 

(9.7) 



NJ 2 



-2 |rp(5) I 2 

5 I 1 AAi I 



+ 4 



Cf + tjt-4 



p-2 \rn(A) I 2 
U .A I 1 AAi I 



cv 



ivj 2 4 
1 



n-2 | T (5) I 2 _| 
u 5 I 1 AAi I ~r 



c F + 



1 \ iV fi 



iv c y 2 i 



p-2 \rpiA) I 2 

I 1 AAi I 



2 \° F + 2N)'f {CsCAyl Re ^ \ W ± 



( 1 



1 



4iV? \4JV? 4 



^-2| T (<S) I 2 i 
u 5 I 1 AAi I 



1 



1 __ / 1 1 

4iV? UiV? + 4 



2|m(- A ) I 2 



U .A I 1 AAi 



W 



qi 



Here, w q = w q a i a , = Wk a £ a . In deriving f)9.7|) we use the following identities: 



{{t a '^t a2 }{t a \t ai }f {T b )< ai 

({t a '\t a2 }[t a2 ,t ai }) ii '(T b ) a ' iai 
({t ai ,t a2 }t b {t a2 ,t ai }) u ' = 2 



l\N c 1 



2 + 2N r 



NJ 2 4 
1 



b\ii' 



1 / 1 

4iV" 2 + UiV? + 4 



Njt b f 
1 



6\ii' 



etc. Since the equilibrium plasmino number density is proportional to the identity, the 
commutator terms on the right-hand side of Eq. ()9.3|) vanish. Therefore if the system is in 
conditions, when the off-equilibrium function Sn~ is perturbatively small, then linearized 
Vlasov-Boltzmann equation (|9.fij) . (|9.7jl in a certain sense is exact. Finally by analogy with 
color current induced by color plasmons (Eq. (9.12) in [T]), we can define the color current 
resulting from the color-plasmino number density jjf\x) = (j'q (x), where now 



a f dq fdn eq (q) 



(2tt)3V du,- 



(w- a -w- a ),j(-\x) = -gt a J- 



a f dq fdn eq (q 



(2vr) 3 V dq 



For closing kinetic equation (|9.6j) (and similar equations for the functions w~ a and 
W^ a it remains only to write out the mean field equation defining a change of mean gauge 
field in a system in a self-consistent manner 



W{x)F» u {x)=fc\x) +j<P(x) + J^(x). 
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Here, the color-plasmon induced current jjpix) is defined by Eq. (9.12) in pQ and j^ xt is 
the external current that plays a part of initial color perturbation. It should be noted 
that an existence of a time-space dependence of background gauge field in the system 
leads, generally speaking, to the replacement of dispersion equations (|2.1j) by 

Re *AZ 1 (q°, q; t, x) = 0, Re *A~ ll (k°, k; t, x) = 0, 

slowly depending on x = (t, x) that in turn results in modification of the dispersion 
relations for plasminos and plasmons: u~ = u;~(q, t,x), uj 1 = u l (k,t,x.). Besides a 
color background gauge field removes an degeneracy of spectrum of plasma excitations, 
when the same dispersion relation corresponds to various oscillations with different color 
indices. The consequence of the last property is the fact that the dispersion relations for 
color plasminos and plasmons acquire a matrix character: uj~ = uj 1 = (u lab ). In 

this case the left-hand side of Eq. ()9.2|) should be supplemented by term of the following 
type: 

1 f du-(q,x) dn n -j 

2 1 <9x ' <9q J ' 

which can be considered as additional 'force' term acting on the color quasiparticles. 
We can assume nevertheless that in the weak-field limit background gauge field weakly 
influences on dispersion properties in the system and for the first approximation one can 
successfully neglected by it. 

10 Conclusion 

In this paper within the framework of the hard thermal loops effective theory we have 
presented the formal scheme for deriving the system of the semiclassical Boltzmann equa- 
tions, which describes the space-time evolution of the soft-quark and soft-gluon number 
densities due to their interactions among themselves in hot non-Abeliam plasma. This 
type of the nonlinear interactions forms in itself an important element in the general 
dynamics of the processes occurring in quark-gluon plasma at least in a weak coupling 
regime. For highly excited plasma states, when typical time of relaxation of the hard 
thermal particle distributions is commensurable with the typical time of relaxation of soft 
oscillations or even significantly exceeds it, along with the kinetic equations for hard ther- 
mal particles it should be used the kinetic equations for soft modes of type ()2.2j) . (j2.3j) . 
However, for completeness of the description the right-hand sides of kinetic equations 
()2.2|) . ()2.3|) should be added by collision terms taking into account the scattering pro- 
cesses of soft plasma excitations on hard thermal particles. For moderate level of plasma 
excitations the processes of this type can even play more important role in the dynamics 
of the system in comparison with the above-mentioned processes of self-interaction of soft 
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modes. In our forthcoming papers a formalism of regular calculation of such collision 
terms will be presented to the full extent. 

Furthermore, we would like to discuss briefly how the approach presented here and 
in our early papers, is agreed with the more traditional calculations from the Feynman 
graphs (like the original calculations of the damping rates for soft quasiparticles within 
the framework of the imaginary-time formalism). 

One of the strongest arguments of plausibility of this approach is the fact that damping 
rate of a plasmino at rest |H] calculated on its basis, exactly coincides with the damping 
rate of the standing plasmino derived by Kobes, Kunstater, and Mak ^2] in the framework 
of the Braaten-Pisarski effective theory. The same can be said about the damping rate 
of a plasmon computed within the framework of our approach [23J. Thus remaining in 
the context of the Blaizot-Iancu equations and deriving from them the effective kinetic 
equations for the soft plasma modes, we are able to reproduce the gauge- invariant damping 
rates of the soft quasiparticles (to the leading order in the coupling), which corresponds 
to the damping rates from the resummed perturbation theory fTJ ^] . 

Notice that in the work [2I|, Braaten and Thoma have pointed to an existence of cor- 
respondence of such a kind in calculation of the soft contribution to energy loss of a heavy 
fermion in hot QED plasma. They have shown that the result of calculation using the rig- 
orous imaginary-time formalism and the result employing the alternative field-theoretical 
calculation (that consists in a simple replacement of photon propagator by effective pho- 
ton propagator) coincide among themselves. At the same time they have noted that "...a 
rigorous justification must await the development of resummation methods for the real- 
time formalism that are as powerful as the resummation methods for the imaginary-time 
formalism! 1 . This in full measure can be related to our approach of deriving the effective 
kinetic equations for the soft modes. 

In recent years significant efforts have been undertaken to construction of the hard 
thermal loop resummation technique in the non-equilibrium field theory |2*5| ESj (see, also 
|27j ) allowing one effectively to work not only with the single-particle propagators, but 
also with three- and four-point functions. The methods and approaches suggested in these 
works enable one at least in principle within the framework of standard methods of real- 
time finite temperature field theory, to construct a calculation scheme of kinetic equations 
for soft quasiparticles in QGP. As a guide idea here, it can be used elegant, thoroughly 
considered approach suggested by Blaizot and Iancu in deriving the Boltzmann equation 
for hard transverse gluons in high-temperature plasma |28j . Their derivation relies on a 
gauge- covariant gradient expansion of the Kadanoff-Baym equations for the gluon two- 
point function. The Boltzmann equation has emerged as the quantum transport equation 
to leading order in g for the gauge- covariant fluctuation 5G of a hard gluon propagator. 
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Besides in the above-mentioned paper, the Kadanoff-Baym equations for the off- 
equilibrium propagator of the soft gluon T> fJiU (X, Y), which are formally identical to those 
for a hard gluon propagator G^(X,Y), are written out. These equations for V^X^Y) 
are used in [2H] only to deduce the relation between the off-equilibrium gauge- covariant 
fluctuation V<{k,x) and the gauge-covariant fluctuation of the leading-order soft gluon 
polarization tensor 8Il < {k ) x). The problem of self-interactions of the soft fields is not con- 
sidered here. However, in principal, there is nothing to forbid the use of these equations 
for research of the soft-field dynamics and construction of the relevant transport equations 
within the framework of the scheme suggested by Blaisot and Iancu. Here, by 511 < (k,x) 
we mean the fluctuation of the next-to-leading order of the soft-gluon self-energy involving 
three- and four-gluon off-equilibrium vertices with soft external lines. The relevant effects 
of self-interaction of the soft fields are encoded in these vertex functions. It would allow us 
to derive the transport equations for the soft gluons directly from the underlying quantum 
field theory and compare them with the equation obtained in this paper in the context 
of the semiclassical approximation. For deriving the kinetic equations for soft fermion 
modes it should be used the Kadanoff-Baym equations for the quark two-point function. 
The development of this approach is also needed to specify the limits of validity of the 
semiclassical kinetic approach suggested in present work to the research of the processes 
of nonlinear interaction of the soft fields in hot QCD plasma. 

Unfortunately on this way considerable technical complications resulting from the dou- 
bling of degrees of freedom arise. Although it is suggested a number of receptions 25 , 
which to a certain extent simplify intermediate computations in deriving convolutions of 
propagators and vertices, nevertheless a practical application of these methods to solving 
specific problem on deriving soft-quasiparticle kinetic equations is at present very com- 
plicated and requires using the symbolic manipulation program. In this sense our less 
strict approach gives more economic method for calculation of desired kinetic equations, 
although it also is very complex and lengthy. The "memory" of initial rigorous theory 
here, is a presence of different types of vertices *T ( - G ^ and *T^' connected with different 
chronological orders of incoming lines. However their amount is small that gives a possi- 
bility to perform calculation in visible form without resorting to a program for evaluating 
real time Feynman amplitudes. It can be assumed that in rigorous approach the final 
kinetic equations (at least their linearized versions) will be practically the same, as the 
kinetic equations obtained within the framework of our more phenomenological approach. 
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Appendix A 



The explicit form of the HTL-induced vertices between quark pair and two gluons is 
defined by the following expressions: 

m(n „i„ , , 2 rdn vav v i> ( t a t b th a \ 

5T$ ab (k 1 ,k 2 ;q 1 ,q 2 ) = -u 2 —- ^ - -. —r r~ r + . —j-. , 

p J 4-7T (v ■ q± — ie)(v ■ q 2 - it) \v ■ (qi + fci) + it v ■ (qi + k 2 ) - it J 

(A.l) 

c-^ro^h/, , x -ifdSl v u v v i> ( t a t h t b t a \ 

5Tf l ) ab {k 1 M-q u q 2 ) = -ul —- — ^ - —. r^— + : ZTTX^ " 

p 7 47r (d • q 1 + ze)(u • q 2 - it) \v ■ (qi + ki) +it v ■ (qi + k 2 ) + it J 

(A.2) 

Below we list the properties of HTL-resummed three-gluon, two-quark -one-gluon and 
two-quark - two-gluon vertex functions, which used in the text 

(*r //WM2 (— k\ — k 2 , ki, k 2 )) = — *r MMlM2 (A;i + k 2 , —k±, —k 2 ) = *V ^^{ki + k 2 , —k 2 , —ki), 

7°rrW(fc;?i,«,))V = *rW(k;q 2 , qi ) = *Tf\-k-- qi -q 2 ), (A.3) 
7 °(Tf (fc;?i,ft)) t 7°= *rf (_fc ; - (Zl ,- fe )= Tf (-fcj-ga,-^), 

(5r«^(fci, A; 2 ; fe ) = arjg^-to, -fc l5 -g 2 , = 5r^ ba (A; 2 , h; Qi, q 2 ), (A.4) 
7°(*r<S )a *(*i, A; 2 ; ?2 )) V = - arjg^fo, A; 2 ; g 2 , ?1 ) = - ST^-h, -k 2 ; -q u -q 2 ), 

ST^fa, k 2 ; q u q 2 ) = - ST ( ^} ba (k 1: k 2 ; q 2 , (A.5) 
-fiST^ih, k 2 ; qi , q 2 ))h° = ~ SF^i-h, -k 2 - -q ± , -q 2 ). 



Appendix B 

In this Appendix we give the explicit expressions for the coefficient functions entering 
into integrands of the conjugate effective currents and sources. For the conjugate effective 
current 

~jf ™*(A<°\ ^<°>) = - gfv^Zf (k, -h; -^^^(^^'("(ft)^^) 

x S(k — k\—q\ — q 2 )dk\dqidq 2 (B.l) 
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we have 

*r m , Q ^ (A:, -fei; -gi, -g 2 ) = STfJy^i-k, fc a ; g x , g 2 ) (B.2) 

- [t% t a f (-*, fc - fci, fc a ) ^(-fc + fci) T^ a (-ft - g 2 ; ft, g 2 ) 

- (t a H a y* T^ 7 (-fc i; -g 2 , + g 2 ) *S lY (-k + ft) T^-fc; k - ft, ft) 

+ (t a t a T k - g 2 , g 2 ) *S 7 y (-A* - ft) T«J ya (-fc a ; -ft, ft + A*). 

Note that this expression cannot be obtained by replacements a # ai, a # /3, —A;, 
fci ^ — fci, and gi j2 ^ — gi. 2 from expression (|4.8jl by virtue of nontrivial (spinor) structure 
of two last terms. However, here we can point to the following formula connecting these 
effective amplitudes 

1 Wl,aP _<?2 ' ~~ ~~ 1 {-KiKli<l2,qi)- 

For Dirac conjugate effective source fy^' 1 ^ ()3.13|) we have expression for conjugate 
effective amplitude in integrand 

•^^(fci.fc; = -{5r£SS 4l< (-fci,-*%i; (B.3) 

+ [t* T<%(g - g i; -g, gi ) - fo) *^ lfl2 (-h ~ h, k u k 2 ) 

- (t a H a ^ T^(A%,; -ft - fc 2 , gi ) (-ft - k 2 ) T«^ a (* i; -g, g - A*) 

- (^)* T^(* i; -ft - A*, ft) *5- 7y (-ft - h) T«^ a (fe; -g, g - A*)}. 

As in the case of effective amplitude (B.2) this effective one cannot be obtained from 
expression (|4.9jl by a simple replacement of indices and momenta of the following type: 
a\ # a 2 , i # ii, A;i # — fci etc. Effective amplitude (B.3) possesses a property 

~(Q)a\a2,ii\ ~{Q)a,2a\,ii\ 

* r M iM2,a/3 (h,k 2 ;qi,-q) = *r M2Ml;Q/3 (k^hiq^-q). (B.4) 

Finally we consider the conjugate effective source fja , Eq. (|3.14|) . By direct calcula- 

tion it is not difficult to show that the conjugate effective amplitude T aaia2a3 is associated 
with r^m^ag by the following expression: 

* T aaia2 «Mi 9i, -92, -gs) = - ^^(-g, -gi, g 2 , g 3 ). (B.5) 
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